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Abstract. We extend the relation between random matrices and 
free probability theory from the level of expectations to the level 
of all correlation functions (which are classical cumulants of traces 
of products of the matrices) . We introduce the notion of "higher 
order freeness" and develop a theory of corresponding free cumu- 
lants. We show that two independent random matrix ensembles 
are free of arbitrary order if one of them is unitarily invariant. We 
prove R-transform formulas for second order freeness. Much of 
the presented theory relies on a detailed study of the properties of 
"partitioned permutations" . 



1. Introduction 

Random matrix models and their large dimension behavior have been 
an important subject of study in Mathematical Physics and Statistics 
since Wishart and Wigner. Global fluctuations of the eigenvalues (that 
is, linear functionals of the eigenvalues) of random matrices have been 
widely investiga ted in the last decade; see, e.g., |Joh98| IDia03t IRad04| 



IMNn4[ lMSSn4j . Roughly speaking, the trend of these investigations 
is that for a wide class of converging random matrix models, the non- 
normalized trace asymptotically behaves like a Gaussian variable whose 
variance only depends on macroscopic parameters such as moments. 
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The philosophy of these results, together with the freeness results of 
Voiculescu served as a motivation for our series of papers on second 
order freeness. 

One of the main achievements of the free probability theory of Voi- 
culescu |Voi9H IVDN92j was an abstract description via the notion of 
"freeness" of the expectation of these Gaussian variables for a large 
class of non-commuting tuples of random mat rices. 

In the previous articles of this series |MS04t IMSS04 we showed that 
for many interesting ensembles of random matrices an analogue of the 
results of Voiculescu for expectations holds also true on the level of 
variances as well; thus pointing in the direction that the structure of 
random matrices and the fine structure of their eigenvalues can be 
studied in much more detail by using the new concept of "second order 
freeness". One of the main obstacles for such a detailed study was 
the absence of an effective machinery for doing concrete calculations 
in this framework. Within free probability theory of first order, such 
a machinery was provided by Voiculescu with the concept of the R- 
transform, and by Speicher with the concept of free cumulants; see, 
e.g., |VDN92[ [N^. 

One of the main achievements of the present article is to develop a 
theory of second order cumulants (and show that the original definition 
of second order freeness from Part I of this series |MS04j is equivalent 
to the vanishing of mixed second order cumulants) and provide the 
corresponding /^-transform machinery. 

In Section 121 we will give a more detailed (but still quite condensed) 
survey of the connection between Voiculescu's free probability theory 
and random matrix theory. We will there also provide the main mo- 
tivation, notions and concepts for our extension of this theory to the 
level of fluctuations (second order), as well as the statement of our 
main results concerning second order cumulants and i?-transforms. 

Having first and second order freeness it is, of course, a natural 
question whether this theory can be generalized to higher orders. It 
turns out that this is the case, most of the general theory is the same for 
all orders. So we will in this paper consider freeness of all orders from 
the very beginning and develop a general theory of higher order freeness 
and higher order cumulants. Let us, however, emphasize that first and 
second order freeness seem to be more important than the higher order 
ones. Actually, we can prove some of the most important results (e.g. 
the i?-transform machinery) only for first and second order, mainly 
because of the complexity of the underlying combinatorial objects. 
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The basic combinatorial notion behind the (usual) free cumulants 
are non-crossing partitions. Basically, passage to higher order free cu- 
mulants corresponds to a change to multi-annular non-crossing per- 
mutations |MJNn4j . or more general objects which we call "partitioned 
permutations" . For much of the conceptual framework there is no dif- 
ference between different levels of freeness, however for many concrete 
questions it seems that increasing the order makes some calculations 
much harder. This relates to the fact that n-th order freeness is de- 
scribed in terms of planar permutations which connect points on n dif- 
ferent circles. Whereas enumeration of all non-crossing permutations 
in the case of one circle is quite easy, the case of two circles gets more 
complicated, but is still feasible; for the case of three or more circles, 
however, the answer does not seem to be of a nice compact form. 

In the present paper we develop the notion and combinatorial ma- 
chinery for freeness of all orders by a careful analysis of the main ex- 
ample: unitarily invariant random matrices. We start with the calcu- 
lation of mixed correlation functions for random matrices and use the 
structure which we observe there as a motivation for our combinatorial 
setup. In this way the concept of partitioned permutations and the 
moment-cumulant relations appear quite canonically. 

We want to point out that even though our notion of second and 
higher order freeness is modeled on the situation found for correlation 
functions of random matrices, this notion and theory also have some 
far-reaching applications. L et us m ention in this respect two points. 



Firstly, recently one of us |Sni05j developed a quite general theory for 



fluctuations of characters and shapes of random Young diagrams con- 
tributing to many natural representations of symmetric groups. The 
results presented there are closely (though, not explicitly) related to 
combinatorics of higher order cumulants. This connection will be stud- 



ied in detail in the part IV of this series [CMSSj where we prove that 
under some mild technical conditions Jucys-Murphy elements, which 
arise naturally in the study of symmetric groups, are examples of free 
random variables of higher order. 

In another direction, the description of subfactors in von Neumann 
algebras via planar algebras |Jon99j relies very much on the notions of 
annular non-crossing partitions and thus resembles the combinatorial 
objects lying at the basis of our theory of second order freeness. This 
indicates that our results could have some relevance for subfactors. 

Overview of the article. In Section |21 we will give a compact survey 
of the connection between Voiculescu's free probability theory and ran- 
dom matrix theory, provide the main motivation, notions and concepts 
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for our extensfon of this theory to the level of fluctuatfons (second or- 
der), as well as the statement of our main results concerning second 
order cumulants and i?-transforms. We will also make a few general 
remarks about higher order freeness. 

In Section 121 we will introduce the basic notions and relevant results 
on permutations, partitions, classical cumulants, Haar unitary random 
matrices, and the Weingarten function. 

In Section m we study the correlation functions (classical cumulants 
of traces) of random matrix models. We will see how those are related 
to cumulants of entries of the matrices for unitarily invariant random 
matrices and we will in particular look on the correlation functions 
for products of two independent ensembles of random matrices, one 
of which is unitarily invariant. The limit of those formulas if the size 
N of the matrices goes to infinity will be the essence of what we are 
going to call "higher order freeness". Also our main combinatorial 
objects, "partitioned permutations", will arise very naturally in these 
calculations. 

In Section El we will forget for a while random variables and just look 
on the combinatorial essence of our formulas, thus dealing with multi- 
plicative functions on partitioned permutations and their convolution. 
The Zeta and Mobius functions on partitioned permutations will play 
an important role in these considerations. 

In Section ini we will derive, for the case of second order, the analogue 
of the R-transform formulas. 

In Section [71 we will finally come back to a (non-commutative) prob- 
abilistic context, give the definition and work out the basic properties 
of "higher order freeness" . 

In Section jSl we introduce the notion of "asymptotic higher order 
freeness" . We show that the Itzykson-Zuber integral encodes all infor- 
mation about higher order freeness. We also indicate how our tech- 
niques can give some insight into the computation of some limits of 
matrix integrals. 

In an appendix. Section El we provide a graphical interpretation of 
partitioned permutations as a special case of "surfaced permutations" . 



2. Motivation and Statement of our Main Results 
Concerning Second Order Freeness and Cumulants 

In this section we will first recall in a quite compact form the main 
connection between Voiculescu's free probability theory and questions 
about random matrices. Then we want to motivate our notion of second 
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order freeness by extending these questions from the level of expecta- 
tions to the level of flu ctuation s. We will recall the relevant results 
from the papers jMS04t IMSS04 and state the main new results of the 
present paper. Even though in the later parts of the paper our treat- 
ment will include freeness of arbitrarily high order, we restrict ourselves 
in this section mainly to the second order. The reason for this is that 
(apart from first order) second order freeness seems to be the most 
important order for applications, so that it seems worthwhile to spell 
out our general results for this case more explicitly. Furthermore, it is 
only there that we have an analogue of i?-transform formulas. We will 
make a few general remarks about higher order freeness at the end of 
this section. 



2.1. Moments of random matrices and asymptotic freeness. 

Assume we know the eigenvalue distribution of two matrices A and B. 
What can we say about the eigenvalue distribution of the sum A + B 
of the matrices? Of course, the latter is not just determined by the 
eigenvalues of A and the eigenvalues of B, but also by the relation 
between the eigenspaces of A and of B. Actually, it is a quite hard 
problem (Horn's conjecture) — which was only solved recently — to 
characterize all possible eigenvalue distributions of A-\- B. However, if 
one is asking this question in the context of A^ x A^-random matrices, 
then in many situations the answer becomes deterministic in the limit 
A^ oo. 

Definition 2.1. Let A = (y4jv)ArGN be a sequence of A^ x A^-random 
matrices. We say that A has a limit eigenvalue distribution if the limit 
of all moments 

an := lim ^[tr(A^)] {n G N) 

exists, where E denotes the expectation and tr the normalized trace. 

In this language, our question becomes: Given two random matrix 
ensembles of A^x A^-random matrices, A = {A^)iy^^ and B = {Biy)^^f^, 
with limit eigenvalue distribution, does also their sum C = {Cj\f)]\f^^, 
with Cjy = An + Bj^, have a limit eigenvalue distribution, and fur- 
thermore, can we calculate the limit moments of C out of the limit 
moments {a^)k>i of A and the limit moments {a^)k>i of -B in a deter- 
ministic way. It turns out that this is the case if the two ensembles are 
in generic position, and then the rule for calculating the limit moments 
of C are given by Voiculescu's concept of "freeness" . Let us recall this 
fundamental result of Voiculescu. 
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Theorem 2.2 (Voiculescu |Voi91p . Let A and B he two random ma- 
trix ensembles of N x N -random matrices, A = (Ajv)7V6n md B = 
Nm> each of them with a limit eigenvalue distribution. Assume 
that A and B are independent (i.e., for each G N, all entries of 
An are independent from all entries of B^), and that at least one of 
them is unitarily invariant (i.e., for each N, the joint distribution of 
the entries does not change if we conjugate the random matrix with an 
arbitrary unitary N x N matrix). Then A and B are asymptotically 
free in the sense of the following definition. 

Definition 2.3 (Voiculescu |Voi85j ). Two random matrix ensembles 
A = {AN)i\i(zfq and B = {BN)Nm with limit eigenvalue distributions are 
asymptotically free if we have for all p > 1 and all n(l), m(l), . . . , n{p), 
m(j)) > 1 that 



lim E 

W— >oo 







One should realize that asymptotic freeness is actually a rule which 
allows to calculate all mixed moments in A and B, i.e. all expressions 

lim ^[tr(A"(i)5"^(^M"(2)5"^(2) . . . 

out of the limit moments of A and the limit moments of -B. In par- 
ticular, this means that all limit moments of A -\- B (which are sums 
of mixed moments) exist, thus A + B has a limit distribution, and are 
actually determined in terms of the limit moments of A and the limit 
moments of B. The actual calculation rule is not directly clear from 
the above definition but a basic result of Voiculescu shows how this can 
be achieved by going over from the moments q;„ to new quantities 
In |Spe94| , the combinatorial structure behind these k„ was revealed 
and the name "free cumulants" was coined for them. Whereas in the 
later parts of this paper we will have to rely crucially on the combina- 
torial description and their extensions to higher orders, as well as on 
the definition of more general "mixed" cumulants, we will here state 
the results in the simplest possible form in terms of generating power 
series, which avoids the use of combinatorial objects. 

Definition 2.4 (Voiculescu |Voi86j . Speicher |Spe94| ). Given the mo- 
ments (an)n>i of some distribution (or limit moments of some random 
matrix ensemble), we define the corresponding free cumulants {Kn)n>i 
by the following relation between their generating power series: If we 
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put 



n>l n>l 



then we require as a relation between these formal power series that 

C{xM{x)) = M{x). 

Voiculescu actually formulated the relation above in a slightly dif- 
ferent way using the so-called i?-transform TZ{x), which is related to 
C{x) by the relation 

C{x) = 1 + zn{x) 

and in terms of the Cauchy transform G{x) corresponding to a measure 
with moments q;„, which is related to M{x) by 

Gix) = —. 

X 

In these terms the equation C{xM{x)) = M{x) says that 

(1) -^ + n{G{x)) = x, 

G{x) 

i.e., that G{x) and K{x) '■= ^ + TZi^x) are inverses of each other under 
composition. 

One should also note that the relation G{xM{x)) = M{x) determines 
the moments uniquely in terms of the cumulants and the other way 
around. The relevance of the Kn and the /^-transform for our problem 
comes from the following result of Voiculescu, which provides, together 
with dH), a very efficient way for calculating eigenvalue distributions of 
the sum of asymptotically free random matrices. 

Theorem 2.5 (Voiculescu |Voi86j ) . Let A and B be two random matrix 
ensembles which are asymptotically free. Denote by k^, , the 
free cumulants of A, B, A + B, respectively. Then one has for alln > 1 
that 

Alternatively, 

7^^+^(a;) =7^^(a;) +7^^(x). 

This theorem is one reason for calling the k„ cumulants, but there 
is also another justification for this, namely they are also the limit of 
classical cumulants of the entries of our random matrix, in the case that 
this is unitarily invariant. This description will follow from our formulas 
(PHj) and (jHUj) . We denote the classical cumulants by kn, considered as 
multi-linear functionals in n arguments. 
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Theorem 2.6. Let A = (Ajv)AreN be a unitarily invariant random, 
matrix ensemble of N x N random matrices Ajsf whose limit eigenvalue 
distribution exists. Then the free cumulants of this matrix ensemble 
can also be expressed as the limit of special classical cumulants of the 
entries of the random matrices: If = {a'l^^)^j^i, then 

_ 1- ]un-l]^ ( (N) (N) (N) . 

for any choice of distinct i(l), . . . , i(n). 

2.2. Fluctuations of random matrices and asymptotic second 
order freeness. There are many more refined questions about the Hm- 
iting eigenvalue distribution of random matrices. In particular, ques- 
tions around fluctuations have received a lot of interest in the last 
decade or so. The main motivation for introducing the concept of "sec- 
ond order freeness" was to understand the global fluctuations of the 
eigenvalues, which means that we look at the probabilistic behavior of 
traces of powers of our matrices. The limiting eigenvalue distribution, 
as considered in the last section, gives us the limit of the average of 
this traces. However, one can make more refined statements about 
their distributions. Consider a random matrix A — {A]sr)Nef>i and look 
on the normalized traces tr(A^). Our assumption of a limit eigenvalue 
distribution means that the limits := limTv^oo -E'[tr(A^)] exist. It 
turned out that in many cases the fluctuation around this limit, 

tr(A^) - ak 

is asymptotically Gaussian of order i.e., the random variable 

N ■ {tr{A%) - ak) = Tr(^^) - Na^ = 'IV(^^ - a^l) 

(where Tr denotes the unnormalized trace) converges for iV — > cxo to 
a normal variable. Actually, the whole family of centered unnormal- 
ized traces {TY{A^) — Nak)k>i converges to a centered Gaussian family. 
(One should note that we restrict all our considerations to complex ran- 
dom matrices; in the case of real random matrices there are additional 
complications, which will be addressed in some future investigations.) 
Thus the main information about fluctuations of our considered en- 
semble is contained in the covariance matrix of the limiting Gaussian 
family, i.e., in the quantities 

am,n:= lim cov(Tr(A^), Tr(A^;^)). 

N^oo 

Let us emphasize that the and the a^.n are actually limits of classical 
cumulants of traces; for the first and second order, with expectation as 
first and variance as second cumulant, this might not be so visible, but 
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it will become evident when we go over to higher orders. Nevertheless, 
the a's will behave and will also be treated like moments; accordingly 
we will call the am,n 'fluctuation moments'. We will later deflne some 
other quantities ,1, which take the role of cumulants in this context. 

This kind of convergence to a Gaussian family was formalized in 
|MSn4j as follows. Note that convergence to Gaussian means that all 
higher order classical cumulants converge to zero. As before, we denote 
the classical cumulants by so ki is just the expectation, and k2 the 
covariance. 

Definition 2.7. Let A = (A7v)ArGN be an ensemble of N x N random 
matrices A^r. We say that it has a second order limit distribution if for 
all m,n >1 the limits 

an := lim /ci(tr(A^)) 

and 

a^,„:= lim fc2(Tr(A^), Tr(A^)) 

N—*oo 

exist and if 

lim fc,(Tr(Af)),...,Tr(Af)))=0 

for all r > 3 and all n{l), . . . , n{r) > 1. 

We can now ask the same kind of question for the limit fluctuations as 
for the limit moments; namely, if we have two random matrix ensembles 
A and B and we know the second order limit distribution of A and the 
second order limit distribution of B, does this imply that we have a 
second order limit distribution for A + B, and, if so, is there an effective 
way for calculating it. Again, we can only hope for a positive solution 
to this if A and B are in a kind of generic position. As it turned out, 
the same requirements as before are sufficient for this. The rule for 
calculating mixed fluctuations constitutes the essence of the deflnition 
of the concept of second order freeness. 

Theorem 2.8 (Mingo, Sniady, Speicher [MSS04j). Let A and B be two 

random matrix ensembles of N x N -random matrices, A = {Am)^^^ 
and B = {Bn)^^^^, each of them having a second order limit distribu- 
tion. Assume that A and B are independent and that at least one of 
them is unitarily invariant. Then A and B are asymptotically free of 
second order in the sense of the following definition. 

Definition 2.9 (Mingo, Speicher |MS04 j). Consider two random ma- 
trix ensembles A = {A^)^^^ and B = {Biy)]y^^, each of them with a 
second order limit distribution. Denote by 

Fjv(n(l), m(l), . . . , n(p), m(p)) 
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the random variable 

The random matrices A = {An)^^^ and B = {B]\r)]\r^^ are asymptoti- 
cally free of second order if for all n,m > 1 

hm h{TT{A% - a^l),TT{B^ - a^l)) = 

and for all p,q > I and . . . , n{p),m{l), . . . , m(p),n(l), . . . , n{q), 

m(l), . . . , m{q) > 1 we have 

lim k2(YN{n{l),m{l), . . . ,n{p),m{p)) ,YN{n{l),m{2), . . . , h{q) , Th{q)) 

N — *oo \ 

ii p ^ q, and otherwise (where we count modulo p for the arguments 
of the indices, i.e., n{i + p) = n{i)) 



lim /c2(YAr(n(l),m(l), . . . ,n{p),m{p)) ,YN(n{p),m{p), . . . ,n(l),m(l) 

P P 

~ ^ ] i^n{i+k)+n{i) ~^n{i+k)'^n{i)) i^m(i+k)+m,{i+l) ~^m{i+k)^m{i+l)) ■ 
k=l i=l 

Again, it is crucial to realize that this definition allows one (albeit in 
a complicated way) to express every second order mixed moment, i.e., 
a limit of the form 

lim A;2(Tr(A^(')E™(') ■ ■ ■ A^^^^E^^^^), Tt{A%^'^bP'^ ■ ■ ■ A'^^'^Bp"^)) 

in terms of the second order limits of A and the second order limits 
of B. In particular, asymptotic freeness of second order also implies 
that the sum A + B of our random matrix ensembles has a second 
order limit distribution and allows one to express them in principle in 
terms of the second order limit distribution of A and the second order 
limit distribution of B. As in the case of first order freeness, it is not 
clear at all how this calculation of the fiuctuations of A + B out of the 
fiuctuations of A and the fluctuations of B can be performed effectively. 
It is one of the main results of the present paper to achieve such an 
effective description. We are able to solve this problem by providing 
a second order cumulant machinery, similar to the first order case. 
Again, the idea is to go over to quantities which behave like cumulants 
in this setting. The actual description of those relies on combinatorial 
objects (annular non-crossing permutations), but as before this can 
be reformulated in terms of formal power series. Let us spell out the 
definition here in this form. (That this is equivalent to our actual 
definition of the cumulants will follow from Theorem 16.31 ) 
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Definition 2.10. Let (a„)„>i and {am,n)m,n>i describe the first and 
second order limit moments of a random matrix ensemble. We define 
the corresponding first and second order free cumulants (k„)„>i and 
{.i^m,,n)m,n>i by the following requirement in terms of the corresponding 
generating power series. Put 

C{x) := 1 + ^KnX'^, C{x,y) := ^ K^,nX"'y'' 

n>l m,n>l 

and 

n>l m,n>l 

Then we require as relations between these formal power series that 

(2) C{xM{x)) = M{x) 
and for the second order 

(3) M(.,.)^^(.MW,,M,.)).i™.il0>l, 
where 

(4) H(..,):^C(..,)-.,4.og(£^M^- 
or equivalently, 

(5) M{x,y) = C{xM{x),yM{y)) ■ ^^^^^"^^^ - 4^^^^^^) 



+ xy 



M{x) M{y) 

dx ^ 



iixMix)) ■ j-iyM{y)) 



{xM{x) -yM{y) f {x-yfJ' 

As in the first order case, instead of the moment power series M(x, y) 
one can consider a kind of second order Cauchy transform, defined by 



G{x,y) :-- 



xy 

If we also define a kind of second order R transform TZ{x, y) by 

T^{.x,y) := —C{x,y), 
xy 

then the formula (jH]) takes on a particularly nice form: 
(6) 
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G{x) is here, as before, the first order Cauchy transform, G{x) = 
lM{l/x). 

The Km^n defined above deserve the name "cumulants" as they hn- 
earize the problem of adding random matrices which are asymptoti- 
cally free of second order. Namely, as will follow from our Theorem 
I7.15[ we have the following theorem, which provides, together with (jH)), 
an effective machinery for calculating the fluctuations of the sum of 
asymptotically free random matrices. 

Theorem 2.11. Let A and B he two random matrix ensembles which 
are asymptotically free. Then one has for all m,n > 1 that 

A+B _ A . B 1 — -I- 

-|- ana i^rn,n ""m,7i ' '^m,n- 

Alternatively, 



n^+^(x) = n^(x) + n''(x) 



and 



7^^+^(x,y) = 7^^(x,y)+7^^(a;,y). 



Again, one can express the second order cumulants as limits of clas- 
sical cumulants of entries of a unitarily invariant matrix. In contrast 
to the first order case, we have now to run over two disjoint cycles in 
the indices of the matrix entries. This theorem will follow from our 
formulas ^ and pijl . 

Theorem 2.12. Let A = {An)n&'n be a unitarily invariant random 
matrix ensemble which has a second order limit distribution. Then 
the second order free cumulants of this matrix ensemble can also be 
expressed as the limit of classical cumulants of the entries of the random 
matrices: If An = {a-if^)fj=i, then 

i^A Km AT^+n-h, 

'^m,n — nniAT^oo "^m+n j(2) ; "i(2)i(3) ' • • • ' "'j(m),i{l)' 

"j(l)i(2)'"i(2)i(3)'---'"i(n),i(l)^ 

for any choice of distinct . . . , i{m),j{l), . . . ,j{n). 

This latter theorem makes it quite obvious that the second order 
cumulants for Gaussian as well as for Wishart matrices vanish iden- 
tically, i.e., 7l{x,y) = and thus we obtain in these cases that the 
second order Cauchy transform is totally determined in terms of the 
first order Cauchy transform (i.e., in terms of the limiting eigenvalue 
distribution) via 

G'ix)G'{y) 1 



(7) Gix,y) 



iG{x) - G{y)f (x - yy 
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This formula for fluctuations of Wishart matrices was also derived by 
Bai and Silverstein in |BS04j . 

2.3. Higher order freeness. The idea for higher order freeness is 
the same as for second order one. For a random matrix ensemble 
A = {Aiy)N(=N we define r-th order limit moments as the scaled limit 
of classical cumulants of r traces of powers of our matrices, 

hm iV2-fc,(Tr(/lf )),..., Tr(/lf))). 

(The choice of A^^~^' is motivated by the fact that this is the leading 
order for many interesting random matrix ensembles, e.g. Gaussian or 
Wishart. Thus our theory of higher order freeness captures the fea- 
tures of random matrix ensembles whose cumulants of traces decay in 
the same way as for Gaussian random matrices.) Then we look at two 
random matrix ensembles A and B which are independent, and one 
of them unitarily invariant. The mixed moments in A and B of order 
r are, in leading order in the limit —>■ oo, determined by the limit 
moments of A up to order r and the limit moments of B up to order 
r. The structure of these formulas motivates directly the definition of 
cumulants of the considered order. The definition of those is in terms 
of a moment-cumulant formula, which gives a moment in terms of cu- 
mulants by summing over special combinatorial objects, which we call 
"partitioned permutations". Most of the theory we develop relies on 
an in depth analysis of properties of these partitioned permutations 
and the corresponding convolution of multiplicative functions on parti- 
tioned permutations. Our definition of "higher order freeness" is then 
in terms of the vanishing of mixed cumulants. It follows quite eas- 
ily that in the first and second order case this gives the same as the 
relations in Definitions 12.31 and 12.91 respectively. For higher orders, 
however, we are not able to find an explicit relation of that type. 

This reflects somehow the observation that our general formulas in 
terms of sums over partitioned permutations are the same for all orders, 
but that evaluating or simplifying these sums (by doing partial summa- 
tions) is beyond our abilities for orders greater than 2. Reformulating 
the combinatorial relation between moments and cumulants in terms 
of generating power series is one prominent example for this. Whereas 
this is quite easy for first order, the complexity of the arguments and 
the solution (given in Definition I2.10j) is much higher for second order, 
and out of reach for higher order. 

One should note that an effective (analytic or symbolic) calculation 
of higher order moments of a sum A+B for A and B free of higher order 
relies usually on the presence of such generating power series formulas. 
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In this sense, we have succeeded in providing an effective machinery 
for deahng with ffuctuations (second order), but we were not able to 
do so for higher order. 

Our results for higher orders are more of a theoretical nature. One of 
the main problems we have to address there is the associativity of the 
notion of higher order freeness. Namely, in order to be an interesting 
concept , our definition that A and B are free of higher order should of 
course imply that any function of A is also free of higher order from 
any function of B. Whereas for first and second order this follows 
quite easily from the equivalent characterization of freeness in terms 
of moments as in Definitions 12.31 and 12.91 the absence of such a char- 
acterization for higher orders makes this a more complicated matter. 
Namely, what we have to see is that the vanishing of mixed cumulants 
in random variables implies also the vanishing of mixed cumulants in 
elements from the generated algebras. This is a quite non-trivial fact 
and requires a careful analysis, see section [7| 



3. Preliminaries 

3.1. Some general notation. For natural numbers m, n G N with 
m < n, we denote by [m, n] the interval of natural numbers between m 
and n, i.e., 

[m, n] := {m, m + 1, m + 2, . . . , n — 1, n}. 

For a matrix A = (ajj)f^-^^, we denote by Tr the unnormalized and by 
tr the normalized trace, 

^ 1 
Tr(A):=^a,„ tr(A) := -Tr(A). 

1=1 

3.2. Permutations. We will denote the set of permutations on n el- 
ements by Sn- We will quite often use the cycle notation for such 
permutations, i.e., vr = Z2, • • • , v) is a cycle which sends ik to ik+i 
{k = 1, . . . , r), where v+i = i\- 

3.2.1. Length function. For a permutation vr G S'„ we denote by t^tt the 
number of cycles of vr and by Ivrl the minimal number of transpositions 
needed to write vr as a product of transpositions. Note that one has 

+ = foi' all vr G Sn- 
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3.2.2. Non-crossing permutations. Let us denote by 7^ G Sn the cycle 

7n = (1,2, ...,n). 
For all TT E Sn one has that 

1"^! + |7n7''~"^| <n — 1. 

If we have equality then we call tt non-crossing. Note that this is 
equivalent to 

#7r + #(7„7r-^) + 

If TT is non-crossing, then so are 7n7r~^ and 7r~^7„; the latter is called 
the (Kreweras) complement of tt. 

We will denote the set of non-crossing permutations in Sn by NC{n). 
Note that such a non-crossing permutation can be identified with a 
non-crossing partition, by forgetting the order on the cycles. There is 
exactly one cyclic order on the blocks of a non-crossing partition which 
makes it into a non-crossing permutation. 

3.2.3. Annular non-crossing permutations. Fix m, n e N and denote 
by 7m,n the product of the two cycles 

lm,n = (1, 2, . . . , m)(m -M, m 2, . . . , m n). 

More generally, we shall denote by 7mi,...,mfc the product of the corre- 
sponding k cycles. 

We call a TT G 5'm+n connected if the pair tt and 7^^^ generates 
a transitive subgroup in Sm+n- A connected permutation tt G Sm+n 
always satisfies 

(8) |7r| |7TO,„7r"^| < m-Fn. 

If TT is connected and if we have equality in that equation then we 
call TT annular non- crossing. Note that if tt is annular non-crossing 
then 7m,ri7r~^ is also annular non-crossing. Again, we call the latter 
the complement of vr. Of course, all the above notations depend on the 
pair (m, n); if we want to emphasize this dependency we will also speak 
about (m, n)-connected permutations and (m, n)-annular non-crossing 
permutations. 

We will denote the set of (m, n)-annular non-crossing permutations 
by SNc{m, n). A cycle of a tt G SNc{m, n) is called a through-cycle if it 
contains points on both cycles. Each tt G Sjsicijn^n) is connected and 
must thus have at least one through-cycle. The subset of S^cij^TT^) 
where all cycles are through-cycles will be denoted by S'^'Q<{m,n). 

Again one can go over from annular non-crossing permutations to 
annular non-crossing partitions by forgetting the cyclic orders on cy- 
cles; however, in the annular case, the relation between non-crossing 
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permutation and non-crossing partition is not one-to-one. Since we will 
not use the language of annular partitions in the present paper, this is 
of no relevance here. 

Annular non-crossing permutations and partitions were introduced 
m jMNOi] : there, many different characterizations — in particular, the 
one (jHl) above in terms of the length function — were given. 

3.3. Partitions. We say that V = {Vi, . . . , Vk} is a partition of a set 
[1, n] if the sets Vi are disjoint and non-empty and their union is equal 
to [1, n]. We call Vi, . . . ,Vk the blocks of partition V. 

If V = {Vi, . . . , 14} and W = {Wi, . . . ,Wi} are partitions of the 
same set, we say that V < W if for every block Vi there exists some 
block Wj such that Vi <^ Wj. For a pair of partitions V, W we denote 
by V V W the smallest partition U such that V <U and W <U. We 
denote by 1„ = {[l,n]} the biggest partition of the set 

If TT G Sn is a permutation, then we can associate to tt in a natural 
way a partition whose blocks consist exactly of the cycles of vr; we will 
denote this partition either by Ott G V{n) or, if the context makes the 
meaning clear, just by tt G V{n). 

For a permutation tt G 5^ we say that a partition V is vr-invariant if 
TT preserves each block of V. This means that 0,^ < V (which we will 
usually write just as vr < V). 

If V = {Vi, . . . , Vk} is a partition of the set [1, n] and if, for 1 < z < /c, 
TTj is a permutation of the set Vi we denote by tti x ■ ■ ■ x vTfc G the 
concatenation of these permutations. We say that tt = tti x ■ ■ ■ x vta; is 
a cycle decomposition if additionally every factor tTj is a cycle. 

3.4. Classical cumulants. Given some classical probability space 

P) we denote by E the expectation with respect to the corresponding 
probability measure. 



and by P) the algebra of random variables for which all mo- 

ments exist. Let us for the following put A := P). 

We extend the linear functional E : ^ — C to a corresponding 
multiplicative functional on all partitions by (V G V{n), ai, . . . , a„ G 



E(a): 




A) 
(9) 



Ev[ai, . . . ,a„] := J]^ E[ai, . . . , 




where we use the notation 



E[ai, . . . ,a„|y] := E(aii ■ ■ -OiJ ^^i V = {ii < ■ ■ ■ < is) E V. 
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Then, for V G V{n), we define the classical cumulants as multi- 
hnear functionals on A by 

(10) kv[ai, . . . ,an\ = ^ Ew[ai, • • • , ■ Mobp(„)(W, V), 

wev(n) 
w<v 

where Mobp(„) denotes the Mobius function on V{n) (see |Rot64p . 
The above definition is, by Mobius inversion on V{n), equivalent to 

E(ai---a„)= ^ k^[ai, . . . ,an]. 

The fcjr are also multiplicative with respect to the blocks of V and thus 
determined by the values of 

kn{ai, . . . ,a„) := ki^[ai, . . . ,a„]. 

Note that we have in particular 

ki{a) = E(a) and ^2(^1, ^2) = E(aia2) — E(ai)E(a2). 

An important property of classical cumulants is the following for- 
mula of Leonov and Shiryaev |LS59j for cumulants with products as 
arguments. 

Let m,n E N and 1 < < i{2) < ■■■ < i{m) = n. Define 
U e V{n) by 

W = {(l,...,z(l)),(^(z) + l,...,z(2)),...,(z(m-l) + l,...,z(m))}. 
Consider now random variables ai, . . . , a„ G ^ and define 

: = ai ■ ■ ■ aj(i) 
M ■ = ai(i)+i ■ ■ ■ ai(2) 

A-m ■ '^•t(m— 1)+1 ■ ■ ■ (^i(m)- 

Then we have 

(11) km{Ai,A2,...,Am)= ^ kv[ai, . . . , ttn]. 

vevin) 

The sum on the right-hand side is running over those partitions of n 
elements which satisfy V V W = 1„, which are, informally speaking, 
those partitions which connect all the arguments of the cumulant km, 
when written in terms of the a^. 

Here is an example for this formula; for A;2 (0102, 0304). In order to 
reduce the number of involved terms we will restrict to the special case 
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where E(aj) = (and thus also ki{ai) = 0) for all i = 1,2,3,4. There 
are three partitions vr G 7^(4) without singletons which satisfy 

7rV{(l,2),(3,4)} = l4, 
namely [a""^ -ar'a^ 



I I 



I \ I 



and thus formula (fTT|) gives in this case 

^2 (010-2, 0304) = /c4(ai, 02, 03, 04) 

+ /c2(oi, 04)^2(02, 03) + ^2(01, 03)^2(02, 04)- 

As a consequence of (fTT|) one has the following important corollary: 
If {ai, . . . , an} and {61, ... , 6„} are independent then 

(12) kyv[aibi, . . . ,anbn] = ^ A;v[ai, . . . , a„] • /cy'i&i, • • • , &„]• 



v,v'ev{n) 

VVV' = W 



3.5. Haar distributed unitary random matrices and the Wein- 
garten function. In the following we will be interested in the asymp- 
totics of special matrix integrals over the group U{N) of unitary N x A^- 
matrices. We always equip the compact group IA{N) with its Haar 
probability measure. A random matrix whose distribution is this mea- 
sure will be called a Haar distributed unitary random matrix. Thus the 
expectation E over this ensemble is given by integrating with respect 
to the Haar measure. 

The expectation of products of entries of Haar distributed unitary 
random matrices can be described in terms of a special function on the 
permutation group. Since such considerations go back to Weingarten 
|Wei78j ■ Collins [Col03j calls this function the Weingarten function 
and denotes it by Wg. We will follow his notation. In the following we 
just recall the relevant informat ion abo ut this Weingarten function, for 



more details we refer to |Col03t ICS04 



We use the following definition of the Weingarten function. For 

-K E Sn and N > n we put 

Wg(A^, Tt) = E[mii ■ ■ ■ UnnUln{l) " " ' Unn{n)], 

where U = {uij)fj^i is an N x N Haar distributed unitary random 
matrix. Sometimes we will suppress the dependency on A^ and just 
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write Wg(7r). This Wg(A^, vr) depends only on the conjugacy class of 
TT. General matrix integrals over the unitary group can be calculated 
as follows: 

(13) E[w,.,. 



This formula for the calculation of moments of the entries of a Haar 
unitary random matrix bears some resemblance to the Wick formula 
for the joint moments of the entries of Gaussian random matrices; thus 
we will call (|13|) the Wick formula for Haar unitary matrices. 

The Weingarten function is quite a complicated object, and its full 
understanding is at the basis of ques tions a round Itzykson-Zuber inte- 
grals. One knows (see, e.g., that the leading order in 



1/A^ is given by |7r| + n and increases in steps of 2. 

3.6. Cumulants of the Weingarten function. We will also need 
some (classical) relative cumulants of the Weingarten function, which 
were introduced in IGol03| §2.3]. As before, let Mob-p(„) be the Mobius 
function on the partially ordered set of partitions of [1,^] ordered by 
inclusion. 

Let us first extend the Weingarten function by multiplicative exten- 
sion, for V > TT, by 

Wg(V,7r) := nWg(vrly), 

v&v 

where 7T\y denotes the restriction of vr to the block V & V (which is 
invariant under vr since vr < V). 

The relative cumulant of the Weingarten function is now, for a < 
V < W, defined by 

(14) CvM^)= J2 M5b(W,W)-Wg(W,a) 

uev(n) 
v<u<w 

Note that, by Mobius inversion, this is, for any cr < V < W, equivalent 
to 

(15) Wg(W,a)= J2 Cv,ui^)- 

v<u<w 

In |Gol03l Cor. 2.9] it was shown that the order of Cv,w(c") is at most 
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4. Correlation functions for random matrices 

4.1. Correlation functions and partitioned permutations. Let 

us consider N x A^-random matrices Bi, . . . ,Bn : O, — > MAr(C). The 
main information we are interested in are the "correlation functions" 
of these matrices, given by classical cumulants of their traces, i.e., 

. . . , S„) := K{Tr{B,), . . . , 

Even though these correlation functions are cumulants, it is more ade- 
quate to consider them as a kind of moments for our random matrices. 
Thus, we will also call them sometimes correlation moments. 

We will also need to consider traces of products which arc best en- 
coded via permutations. Thus, for n e Sn, 99 (tt) , . . . , i?„] shall 
mean that we take cumulants of traces of products along the cycles of 
TT. For an n-tuple B = {Bi, . . . , Bn) of random matrices and a cycle 
c = {ii, i2, ■ ■ ■ , ik) with k < n we denote 

B\c :— Bi^Bi^ ' - ' ^ik- 

(We do not distinguish between products which differ by a cyclic rota- 
tion of the factors; however, in order to make this definition well-defined 
we could normalize our cycle c = (ii, ^2, • • • , 'i-k) by the requirement that 
ii is the smallest among the appearing numbers.) For any n G S{n) 
and any n-tuple B = {Bi, . . . , Bn) of random matrices we put 

ip{7r)[Bi, . . . , Bn] := (pr{B\c„ . . . , B\c,), 

where n consists of the cycles Ci, . . . , c^. 
Example: 

ip{{l, 3)(2, 5, 4))[Bi, B2, ^3, ^4, ^5] = ^2iBA B^B.B,) 

^k2{Tr{B^B:,),Tr{B2B,B,)) 

Furthermore, we also need to consider more general products of such 
</9(7r)'s. In order to index such products we will use pairs (V, vr) where 
TT is, as above, an element in Sn, and V G V{n) is a partition which is 
compatible with the cycle structure of tt, i.e., each block of V is fixed 
under tt, or to put it another way, V > tt. In the latter inequality we 
use the convention that we identify a permutation with the partition 
corresponding to its cycles if this identification is obvious from the 
structure of the formula; we will write this partition 0^ or just if no 
confusion will result. 

Notation 4.1. A partitioned permutation is a pair (V, tt) consisting of 
n e Sn and V G V{ri) with V > tt. We will denote the set of partitioned 
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permutations of n elements by VS (n) . We will also put 

VS := [j VS{n). 

neN 

For such a (V, vr) G VS we denote finally 

y.(V, vr)[Si, . . . , S„] := J] ¥^(vr|v)[5i, . . . , 

vgv 

Example: 

</.({l, 3, 4}{2}, (1, 3)(2)(4)) 52, i^a, B^] 

= fc2(Tr(5i53),Tr(54)) ■ A;i(Tr(52)) 

Let us denote by Tr^- as usual a product of traces along the cycles of 
cr. Then we have the relation 

E{Tr^[Ai,...,A„]}= 

wev(n) 

W><T 

By using the formula (lllj) of Leonov and Shiryaev one sees that in 
terms of the entries of our matrices Bk = {h^^^)f.j^^ our (p{U,'y) can 
also be written as 
(17) 

N 

^{U,j)[Bu...,Bn]= J2 ^^\''mih{i)V ■ ■ ■ ^^^Z)ih{n))]- 

v^^u^ .(1),...,W=1 

4.2. Moments of unitarily invariant random matrices. For uni- 
tarily invariant random matrices there exists a definite relation between 
cumulants of traces and cumulants of entries. We want to work out this 
connection in this section. 

Definition 4.2. Random matrices Ai, . . . ,An are called unitarily in- 
variant if the joint distribution of all their entries does not change by 
global conjugation with any unitary matrix, i.e., if, for any unitary ma- 
trix U, the matrix- valued random variables Ai, . . . , An : Q ^ Mjv(C) 
have the same joint distribution as the matrix-valued random variables 
UAiU*, . . . , UAnU* : n Mn{C). 

Let Ai, . . . , An he unitarily invariant random matrices. We will now 
try expressing the microscopic quantities "cumulants of entries of the 
A" in terms of the macroscopic quantities "cumulants of traces of 
products of the A^" . 
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In order to make this connection we have to use the unitary invari- 
ance of our ensemble. By definition, this means that Ai, . . . ,An has the 
same distribution as Ai, . . . , A„ where Ai :— UAiU*. Since this holds 

for any unitary U, the same is true after averaging over such U, i.e., we 
can take in the definition of the Ai the U as Haar distributed unitary 
random matrices, independent from Ai, . . . , An- This reduces calcula- 
tions for unitarily invariant ensembles essentially to properties of Haar 
unitary random matrices; in particular, the Wick formula for the f/'s 
implies that we have an analogous Wick formula for joint moments in 
the entries of the A^. Let us write Ak = {a\j^)f^j^-^ and Ak = {a\j^)f^j^i. 
Then we can calculate: 

r (1) (") \ 

~ ^{'^pih'^rijl ■ ■ ■ '^Pnin'^rnjn} ' ^{^iljl ' ' ' ^injn^ 

= E E W.>.Wg(a7r-i)-E{ai;;---a2j 

i,j TT,a€Sn 

neS„ 

where 
(18) 

^(7r)[A,, . . . , A„] : = 5] Wg(a7r-) • E{a;i^^^ . . . att^^,} 

aeSn i 

= J2 Wg((77r-i) • E{Tr,[A^, A„]}. 

creSn 

= J] Wg((77r-^) • v{yV,a)[A,,...,Ar,] 

aeSn wer(n) 

>V>CT 

= J2 Wg(a7r-^)-^(W,(j)[A,...,A„]. 

{W,a)erS{n) 

The important point here is that Q{7r)[Ai, . . . , An] depends only on the 

macroscopic correlation moments of A. 

We can extend the above to products of expectations by 



veSn 
■7r<V 
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where Q{V,it) is given by multiplicative extension: 
g{V, 7r)[^i, . . . , A] : = n ^(^1^) t^i' • • • ' ^"1^] 



(19) = J] Wg(V,a7r-i)-9^(>V,a)[Ai,...,A„] 



(yV,<T)6PS(n) 
W<V 



Now we can look on the cumulants of the entries of our unitarily 
invariant matrices Ai] they are given by 

Wl, . . . , 4±} = J2 Mobp(„)(W, V) • E^[a(i^) ^, . . . 4±] 



uer(n) 



J2 Yl ^-'P°^ ■ Mobp(„)(W, V) • g{U, 7r)[Ai, ...,An] 

J2^r,po. Yl ^0hrin){l{,V)-g{U,7T)[A,,...,A^]. 



TreS„ UeT(n) 
7r<V V>U>TV 



With the definition 

(20) K{V,7r)[Ai,...,An]:^ J] Mobp(„)(W, V)-g?(W, 7r)[^i, . . . , 



uer(n) 

V>U>TV 



we thereby get 

(21) /cvjaW ^, . . . , ai:lj = 5,,,.^ ■ «:(V, ti)[A,, . . . , 



7rSS„ 
7r<V 



It follows that 

AT 

(^(W,7)[A,...,^n] = Y ^v[4J),(^(i)),...,al"4(7{n))] 

N 



^ Y Y Y ^^°7,i°T • /«(v, tt) [^1, . . . , 

V<W i(l),...,j(„) = l TTGSn 
VV7=W •ir<V 



VV7=W 7r<V 



Since V V 7 = W is, under the assumption tt < V, equivalent to V V 

77r~"'^ = W we can write this also as 

(22) (^(W,7)[A,...,^n]- E '^(V,7r)[A,...,A„]-7V#(7--^). 

(V,w)eP5(n) 
VV77r-l=W 
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Remark 4.3. 1) Note that although the quantity n is defined by (j2(Jj) 
in terms of the macroscopic moments of the Ai, they have also a very 
concrete meaning in terms of cumulants of entries of the A^. Namely, 
if we choose tt ^ Sn and distinct 1 < . . . ,i{n) < N then equation 
becomes, when we set V = 1„, 

(23) n)[Ai, . . . ,An] = ^n(0'l(i)j(7r(l))5 • • • 5 (^l(l)i(TT{n))) 

as the the only term in the sum that survives is the one for vr. 

2) Equation should be considered as a kind of moment- cumulant 
formula in our context, thus it should contain all information for defin- 
ing the "cumulants" k, in terms of the moments (f. Actually, we can 
solve this linear system of equations for k in terms of ip, by using equa- 
tion (pUj) to define k and equation (fT^ for Q. 

K{V,7T)[Ai,...,An] 

= J2 Mobp(„)(W,V)- Yl Wg(W,cr7r-i)-y.(>V,a)[Ai,...,A„] 

V>U>n W<U 

V>M>7rVW 

Thus, by using the relative cumulants of the Weingarten function from 

(HH), we get finally 

(24) 

K{V,n)[Ai,...,An]= Y '^(>V,(t)[Ai,...,A„] ■avw,v(^xvr-i). 

w<v 

3) One should also note that we have defined the Weingarten function 
only for N > n; thus in the above formulas we should always consider 
sufficiently large N. This is consistent with the observation that the 
system of equations ()22|) might not be invertible for too small; the 
matrix (A^#(°"'^ )) „ is invertible for N > n, however, in general 
not for all < n (e.g, clearly not for A^ = 1). One can make sense of 
some f ormulas involving the Weingarten function also for N < n (see 



CS04j ). However, since we are mainly interested in the asymptotic 



behavior of our formulas for A^ — > oo, we will not elaborate on this. 

4.3. Product of two independent ensembles. Let us now calculate 
the correlation functions for a product of two independent ensembles 
Ai, . . . , An and Bi, . . . , Bn of random matrices, where we assume that 
one of them, let's say the -Bj's, is unitarily invariant. We have, by 
using (fTTj) and the special version (|T^ of the formula of Leonov and 
Shiryaev, the following: 
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^(U,-f)[AiBi,...,AnBr. 

.(1) l.t.._.n 

i(l)i(7(l))' • • • ' "i(n)i(7(n))J 



., , VVV'V7=W 



i(l)i(7(l))' • • • ' "i(n)i(7W)J 

^ ^ «:(V,7r)[Ai,...,A„]- 



ij' V,V'<M T&Sn 



W>V>7r 



V'<M J 
V'VVV7=M 



In order to evaluate the second factor we note first that, under the 
assumption vr < V, the condition V V V V 7 = W is equivalent to 

V V V V 7r^^7 = U. Next, we rewrite the sum over all V G V{n) with 

V < W and V'V VV 71-^7 = W as a double sum over all W G V{n) 
with VV W = W and all V G P(n) with V <W and V V n-^^t = W. 



^V'[^i{l)j(7r-l7{l))' • • • ' ^i{n)j(7r-i7(n))] 

'ev(n) i 

V'VVV7=W 

wev(n) v'<w i 

VVW=M v'V7r-l7 = W 

^ y;(W,7r-i7)[i?i,...,5„]. 



yve7'{n) 

W>7r-l7,VVW=W 



Thus we finally get 
<^(W,7)[v4iEi,...,A„5„] 



EE E «:(V,7r)[Ai,...,A„]-y.(>V,7r-S)[5i,...,5n] 

tGS„ V6P{n) wev(,,i) 

U>V>TT w>7r~l7,VVW=M 

5^ «(V, vr) [Ai, . . . , A„] ■ <^(W, TT-S) . . . , S„] . 



{V,7r),(W,<T)6-p5(n) 
VVW=W,7rCT=7 



Let us summarize the result of our calculations in the following the- 
orem. In order to indicate that our main formulas are valid for any 
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fixed N, we will decorate the relevant quantities with a superscript 
Note that up to now we have not made any asymptotic consideration. 

Theorem 4.4. Let Mn '■= (g) L°°{Q) be an ensemble of N x N- 
random matrices. Define correlation functions (pl!^^ on M. ^ by (n & N, 
Di,...,Dne Mn) 

(25) ...,D^):= Kmn,), Tr(D„)) 

and corresponding "cumulant functions" k^^^ (for n < N) by 
(26) 

w<v 

or equivalently by the implicit system of equations 

(27) v^(^)(W, 7) [D,, ...,D^]=Y, k(^)(V, n)[D,, . . . , D„] ■ N*'^^^-'\ 

V,7r 

where the sum is over all V G V{n) all n E Sn such that tt < V and 
VV77r-i = U. 

1) Let An be an algebra of unitarily invariant random matrices in 
M.N- Then we have for all n < N , all distinct i(l), . . . ,i{n), all A^ = 
(o,-?)^- 1 € A, and all vr G S'„ that 

(28) K^^^ln, 7r)[Ai, . . . ,An\ = • • • , «'!(")i('r(n)) ) 

2) Assume that we have two subalgebras An and Bn of M.n such 
that 

o An is a unitarily invariant ensemble, 
o An o-nd Bn are independent. 
Then we have for all n E N with n < N and all Ai, . . . , An G An 
and Bi, . . . ,B„ G Bm ■' 

(29) </.(^)(^^,7)[^i5i,...,A„B„] 

= '^^^HV,vr)[^,...,A„]-y,W(W,a)[5i,...,i?„]. 

V,7r,W,<T 

where the sum is over all V, W G V{n) and all tt, a G 5*^1 such that 
71 <V, a <W, VyW = U, and-f = na. 

4.4. Large N asymptotics for moments and cumulants. Our 

main interest in this paper will be the large N limit of formula ()29|). 
This structure in leading order between independent ensembles of ran- 
dom matrices which are randomly rotated against each other will be 
captured in our abstract notion of higher order freeness. 
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Of course, now we must make an assumption about the asymptotic 
behavior in N of our correlation functions. We will require that the 
cumulants of traces of our random matrices decays in N with the same 
order as in the case of Gaussian or Wishart random matrices. In these 
cases one has very detailed "genus expansions" for those cumulants; 
see, e.g. jOkoOOt IMN04j and one knows that the n-th cumulant of 
unnormalized traces in polynomials of those random matrices decays 
like A^2-n (^ggg g_g_ jMS()4[ Thm. 3.1 and Thm. 3.5]). 

Definition 4.5. Let, for each iV G N, Bf\ B^^'^ C Mj^®L°^-{n) 
he N X A^-random matrices. Suppose that the leading term of the 
correlation moments of b['^\ . . . , Br'^^ are of order 2 — n, i.e., that for 
all n G N and all polynomials pi, . . . ,pt in r non-commuting variables 
the limits 

hm , • • • , Bi^\ . . .MBf \ . . . , 5f ))) ■ AT"- 

exist. Then we will say that {b['^\ . . . , Br'^^} has limit distributions of 
all orders. Let B be the free algebra generated by generators bi, . . . ,br- 
Then we define the limit correlation functions of B by 

= hm • • • , 5f • • ■MBi^'K . . . , 5f ))) ■ iV"- 

N—i-oo 

Note that this assumption implies that the leading term for the 
quantities (f^^\V,7r) is of order 2#(V) — #(vr). Indeed, if V has k 
blocks and the z*'^ block of V contains rj cycles of n then (f^'^\V, tt) = 
■ ■ ■ V^T-fc and each Lp^ has order 2 — r^. Then the order of (f^'^\V, tt) 
is (2-ri) + --- + (2-rfc) = 2A;-(ri + --- + rfc) = 2#(V)-#(7r). Thus 

V9(V,7r)(pi(6i, . . .,br), . . . ,Ptibi, . . . ,6^)) 

= hm ^(^)(V, n){p,iBf\ . . . , . . . • • • , Bi^^)) 

. ^-2#(V)+#(7r) 

From formula ()27|) one can deduce that the leading order of k^^^ (V, tt) 
is given by the term (U, 7) = (V, tt) and thus must be of order 



N- 



-n+2#V-#7r 



(Indeed, this also follows from equation and the leading order of 
the relative cumulant of the Weingarten function given in equation 
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Thus we can define the limiting cumulant functions to be the hmit 
of the leading order of the cumulants by the equation 
(30) 

When (V, vr) = 7^) and Bi = B2 = ■ ■ ■ = = B equation (1^ 
becomes 

Thus to prove Theorem 12 .61 we must show that 7„)[i?, . . . ,B] 

■ N"^"^ converges to the n*^ free cumulant of the hmiting eigenvalue 
distribution of B^^\ 

When (V, vr) = (lm+n,7m,n) equation becomes 

n^m+n, 7m,n)[-B, ■ ■ ■ , B] = ^m+n j(2) 5 • • • 5^i(r!)i(l)) 

Thus to prove Theorem 12.121 we must show that K^'^\lm+n,lm,n)[B, 
. . . ,B] ■ A^™-+" converges to „ the (m, n)*^ free cumulant of second 
order of the limiting second order distribution of B^^\ 

4.5. Length functions. We want to understand the asymptotic be- 
havior of formula ((221) • The leading order in of the right hand side 
is given by 

-n + 2#V - #7r + 2#W - #a = n + (|7r| - 2|V|) + (|(t| - 2|W|), 

whereas the leading order of the left hand side is given by 

2#W - #7 = 2#(V V W) - #((T7r) = n + {\7ra\ - 2|V V 

This suggests the introducing of the following "length functions" for 
permutations, partitions, and partitioned permutations. 

Notation 4.6. 

(1) For V G V{n) and tt G we put 

|V| :=n-#V 
|7r| := n — t^tt. 

(2) For any (V, vr) G VS{n) we put 

|(V,7r)|:=2|V|-|7r|=n-(2#V-#7r). 

Let us first observe that these quantities behave actually like a 
length. It is clear from the definition that they are always non-negative; 
that they also obey a triangle inequality is the content of the next 
lemma. 



Lemma 4.7. 
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(1) For all 7i,a E Sn we have 

|7rcr| < |7r| + |cr|. 

(2) For all V, W G V{n) we have 

|vv wi < |v| + 

(3) For all partitioned permutations (V, vr), (W, a) G VS{n) we 
have 

|(VVW,7ra)| < |(V,7r)| + |(W,a)|. 

Proof. (1) This is well-known, since |7r| is the minimal number of factors 
needed to write tt as a product of transpositions. 

(2) Each block i? of W can glue at most 4i^B — 1 many blocks of V 
together, i.e., W can glue at most n — #W many blocks of V together, 
thus the difference between |V| and |V V >V| cannot exceed n — #>V 
and hence 

#V - #(V V W) < n - #>V. 

This is equivalent to our assertion. 

(3) We prove this, for fixed vr and a by induction over |V| + |W|. The 
smallest possible value of the latter appears for |V| = |7r| and |VV| = |o"| 
(i.e., V = 0^ and W = 0^). But then we have (since VV W > ttct) 

2|V VW| - |7ra| < |VV W| < |V| + 

which is exactly our assertion for this case. For the induction step, on 
the other side, one only has to observe that if one increases |V| (or |W|) 
by one then |V V W| can also increase by at most 1. □ 

Remark 4.8. 1) Note that the triangle inequality for partitioned per- 
mutations together with (j^^ implies the following. Given random ma- 
trices A = (v4Ar)jvgN and B = {B]y)Nm which have limit distributions 
of all orders. If A and B are independent and at least one of them 
is unitarily invariant, then C = {Ciy)Nm with Cjy := AjyB^^ also has 
limit distributions of all orders. 

2) Since we know that Gaussian and Wishart random matrices have 
limit distributions of all orders (see e.g. |MS04| Thm. 3.1 and Thm. 
3.5]), and since they are unitarily invariant, it follows by induction 
from the previous part that any polynomial in independent Gaussian 
and Wishart matrices has limit distributions of all orders. 

4.6. Multiplication of partitioned permutations. Suppose 

. . . , -BI^'' } has limit distributions of all orders. Then the left hand 

side of equation ^I7\ has order A^2#(w)-#(7) g^j^^ ^j^g right hand side of 
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equation ()27|) has order A^~"+2#(v)-#(7r)+|77r Thus the only terms of 
the right hand side that have order A^2#{w)-#(7) ^^^.^ those for which 

2#(W) - #(7) = -n + 2#(V) - #(7r) + \^n-'\ 
i.e. for which |(W,7)| = |(V, 7r)| + |77r^-^|. Hence 

V.(^)(W,7)[5f\...,5f)] 

(V,7r)e-P5(n) 
VV77r-l=W 
|(W,7)NI(V,7r)t + |7,r-l| 

+ o(Ar2#(")-#(7)-2) 

Thus after taking hmits we have 

(31) </.(W,7)[&i,...,M = Yl «^(V,vr)[6i,...,6„] 

(V,7r)e-P5(n) 

where the sum is over all (V, vr) in VS{n) such that V V 77r^^ = U and 
|(W,7)| = |(V,7r)| + |77r-i|. 

A similar analysis of equation (|^^ gives that for independent {A^^\ 
. . . , A\^''} and {B'^\ . . . , -Bn^^} with the ^^-^'''s unitarily invariant and 
both having limit distributions of all orders we have 

^(^)(w,7)[Ar5r\...,4^)5n 

= Y '^^"^HV,vr)[4^\...,AW].^W(VV,a)[5r,...,i?f)] 

(V,7r),(W,CT)S-p5(n) 
|(V,7r)l + |(W,CT)| = |(VVW,ir(T)| 

+ 0(A^2#(W)-#(7)-2) 

and again after taking limits 

(32) v9(W,7)[ai6i, . . . ,a„6„] 

= Y K(V,7r)[ai, . . . ,a„] ■ v5(>V,(t)[6i, . . . , 6„,] 

(V,7r), (W,o-)e-P<S(n) 

where the sum is over all (V, vr), (W, a) G VS{n) such that 

o V V W = W 
o vrcr = 7 

o |(V,7r)| + |(>V,a)| = |(W,7)| 
In order to write this in a more compact form it is convenient to 
define a multiplication for partitioned permutations (in CVSiji)) as 
follows. 
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Definition 4.9. For (V, 7r),(>V, a) G VS{n) we define tlieir product 
as follows. 

(33) (V,7r)-(W,a) : = 

_ f(VV W.vra) if |(V,7r)| + \{W,a)\ = \{VWW,7ia)\, 
1 otherwise. 

Proposition 4.10. The multiplication defined in Definition ^ ■ 9\ is as- 
sociative. 

Proof. We have to check that 

(34) ((V, tt) ■ (W, a)) ■ {U, r) = (V, n) ■ ((W, a) ■ {U, r)) . 

Since both sides are equal to (V V W V W, vrar) in case they do not 
vanish, we have to see that the conditions for non-vanishing are for 
both sides the same. 

The conditions for the left hand side are 

|(V,7r)| + |(W,a)| = |(VVW,7ra)| 

and 

|(VV W,7ra)| + |(W,r)| = |(VV W V W, 7rar)|. 
These imply 

|(V,7r)| + |(W,a)| + \{U,t)\ = |(W V W VW,7rar)| 

< |(V,7r)| + |(>VVW,ar)|, 

However, the triangle inequality 

|(WVW,ar)| < \{W,a)\ + \{U,T)\ 

yields that we have actually equality in the above inequality, thus lead- 
ing to 

|(W,a)| + |(W,r)| = |(WVW,(7r)| 

and 

|(V,7r)| + |(WVW,ar)| = |(VV W V W, Trar)]. 

These are exactly the two conditions for the vanishing of the right hand 
side of dnH). The other direction goes analogously □ 

Now we can write formulas ()31|) and (jH^ in convolution form 

(35) ^(W,7)[6i,...,M = Yl '^(V,7r)[6i,...,6„] 

(V,7r)ePS(n) 
(V,7r)-(0,77r-l) = (M,7) 
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and 

(36) ip(U,'j)[aibi,...,anbn] 

= ^ K{V,7i)[ai, . . . ,a„] ■ (p{yV,a)[bi, . . . , 6„] 

(V,7r)-(W,CT) = {W,7) 

Note that both <f{V, tt) and /t(V, vr) are multiphcative in the sense that 
they factor according to the decomposition of V into blocks. 

The philosophy for our definition of higher order freeness will be 
that equation (j35p is the analogue of the moment-cumulant formula 
and shall be used to define the quantities k, which will thus take on 
the role of cumulants in our theory - whereas the are the moments 
(see Definition 17. 4p . We shall define higher order freeness by requiring 
the vanishing of mixed cumulants, see Definition 17.61 On the other 
hand, equation (jHBj) would be another way of expressing the fact that 
the a's are free from the 6's. Of course, we will have to prove that those 
two possibilities are actually equivalent (see Theorem 17. 9|) . 

5. Multiplicative functions on partitioned permutations 

and their convolution 

5.1. Convolution of multiplicative functions. Formulas ()35|) and 
above are a generalization of the formulas describing first order 
freeness in terms of cumulants and convolution of multiplicative func- 
tions on non-crossing partitions. Since the dependence on the random 
matrices is irrelevant for this structure we will free ourselves in this 
section from the random matrices and look on the combinatorial heart 
of the observed formulas. In Section [7[ we will return to the more gen- 
eral situation involving multiplicative functions which depend also on 
random matrices or more generally elements from an algebra. 

Definition 5.1. 

(1) We denote by VS the set of partitioned permutations on an 
arbitrary number of elements, i.e., 

VS=[j VS{n). 

neN 

(2) For two functions 

we define their convolution 



f*g:VS 



C 
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by 

(/*^)(W,7):= Yl /(V,7r)^(W,a) 

for any (^,7) G VS{n). 

Definition 5.2. A function / : VS — * C is called multiplicative if 
/(l„,7r) depends only on the conjugacy class of tt and we have 

f{V,T^)=\{f{lvMv) 

v&v 

Our main interest will be in multiplicative functions. It is easy to 
see that the convolution of two multiplicative functions is again mul- 
tiplicative. It is clear that a multiplicative function is determined by 
the values of /(!„, tt) for all n e N and all 

An important example of a multiplicative function is the 5-function 
presented below. 

Notation 5.3. The 5-function on VS is the multiplicative function 
determined by 

\, ifn = l 
0, otherwise. 



5(ln,7r) 



Thus for (W, tt) G VS{n) 

1, if (W,7r) = (On, (1)(2) . . . (n)) for some n 
0, otherwise. 



5(W,7r) 



Proposition 5.4. The convolution of multiplicative functions on VS 
is commutative and 5 is the unit element. 

Proof. It is clear that S is the unit element. For commutativity, we 
note that for multiplicative functions we have 

/(V,7r) = /(V,7r-i), 

and thus 

{g * 7) = * 7-') - E 9{V: 7r)/(W, a). 

(V,7r),(VV,a)eP5(n) 
(V,7r)-(W,o-)=(W,7-l) 

Since the condition (V, tt) ■ (W, a) ~ {U, 7"^) is equivalent to the con- 
dition (W, (7~^) • (V, 7r~^) = {U, 7) we can continue with 

{9 * /) (W, 7) = Yl /(W, a-^)^(V, TT-^) ^{f*g) (W, 7) . 

(V,7r),(W,cr)SPS(n) 
(W,CT-l)-(V,7r-l) = (W,7) 

□ 



34 B. COLLINS, J. A. MINGO, P. SNIADY, AND R. SPEICHER 

5.2. Factorizations. Let us now try to characterize the non-trivial 
factorizations (W, 7) = (V, tt) • (W, a) appearing in the definition of our 
convolution. Let us first observe some simple general inequalities. 

Lemma 5.5. 

(1) For permutations tt, cr e S{n) we have 

\tt\ + \cr\ + Itto"! > 2|7rV(T|. 

(2) For partitions V2 < Vi and W2 < Wi we have 

+ |Vi| + IV2 VW2I > |Vi VWil + IW2I + IV2I 

and 

|Vi V W2I + IV2 V Wil > |Vi V Wil + IV2 V W2I. 

Proof. (1) By the triangle inequality for partitioned permutations we 
have 

|(0,V0„7ra)| < |(0,,7r)| + |(0„a)|, 

i.e., 

(37) 2|7r V(t| - \7ra\ < |7r| + \a\. 

(2) Consider first the special case Wi = W2 = W. Then we clearly 
have 

#(V2 V W) - #(Vi V W) < #V2 - #Vi, 

which leads to 

|Vi vw| - IV2 vw| < |Vi| - IV2I. 

From this the general case follows by 

|Vi V Wil - IV2 V W2I = |Vi V Wil - |Vi V W2I + |Vi V W2I - IV2 V W2 

< iWil - IW2I + |Vi| - IV2I. 

The second inequality follows from this as follows: 

|Vi V Wii - |Vi V W2I = |Vi V (V2 V - |Vi V (V2 V >V2)| 

< IV2 V Wi| - IV2 VW2I. 

□ 

Theorem 5.6. For (V,7r), (W, cr) e VS{n) the equation 
(V, tt) • (W, (j) = ( VV W, na) 
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is equivalent to the conjunction of the following four conditions: 
\7t\ + \a\ + \7ia\ = 2\7i V cr|, 
|V| + \n\/ a\ = |7r| + |VV (r\, 
|W| + Itt Va| = |a| + |7r V >V|, 
I VV cr| + |7r VWI = |V VWI + Itt V cr|. 

Proof. Adding the four inequalities given by Lemma f5. 51 

1^1 + 1^1 + K^l — 2|7l' V cr|, 
2|V| + 2|7r Va| > 2|7r| +2|VVa|, 
2|W| +2|7r Va| > 2|(t| + 2|7r V W|, 
2|V Va| + 2|7r V Wl > 2|V V >V| + 2|7r Va| 

gives 

2|V| - |7r| + 2\W\ - \a\ > 2|V V W| - Ivral, 

i.e., 

|(V,7r)| + |(W,a)| > |(VVW,7ra). 

Since (V, tt) • (W, a) = (V V W, vrcr) means that we require equahty in 
the last inequality, this is equivalent to having equality in all the four 
inequahties. □ 

The conditions describing our factorizations have a quite geometrical 
meaning. Let us elaborate on this in the following. 

Definition 5.7. Let 7 G S{n) be a fixed permutation. 

(1) A permutation vr G S{n) is called 'y -planar if 

|7r| + |7r^"^7| + I7I = 2|7rV7|. 

(2) A partitioned permutation (V, tt) G VS{n) is called '-^-minimal 
if 

IVV7I — |7rV7| = |V| — |7r|. 

Remark 5.8. z) It is easy to check (for example, by calculating the Eu- 
ler characteristic) that 7-planarity of tt corresponds indeed to a planar 
diagram, i.e. one can draw a planar graph representing permutations 
7 and TT without any crossings. The most important cases are when 
7 consists of a single cycle |Bia97j and when 7 consists of two cycles 
|MNn4| . 

ii) The notion of 7-minimality of (V, vr) means that V connects only 
blocks of TT which are not already connected by 7. 

Hi) If (V, vr) satisfies both (1) and (2) of Definition IKTIthen (V, 7r)(0, -R-^-i) = 
(1,7), by Theorem 15.61 
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Corollary 5.9. Assume that we have the equation 

(W,7) = (V,7r)-(W,a). 

Then n and a must be ^-planar and (V, tt) and (W, cr) must be 7- 

minimal. 

5.3. Factorizations of disc and tunnel permutations. 

Notation 5.10. i ) We call (V, vr) G VSn a disc permutation if V = 0^; 
the latter is equivalent to the condition |V| = \a\. For tt G Sn, by (0, tt) 
we will always mean the disc permutation 

(0,7r) := (0,,7r) G VS{n). 

a) We call (V, tt) G VSn a tunnel permutation if |V| = |7r| + 1. This 
means that V is obtained from tt by joining a pair of cycles; i.e. one 
block of V contains exactly two cycles of tt and all other blocks contain 
only one cycles of tt. 

A motivation for those names comes from the identification between 
partitioned permutations and so-called surfaced permutations; see the 
Appendix for more information on this. 

Our goal is now to understand more explicitly the factorizations of 
disc and tunnel permutations. (It will turn out that those are the 
relevant ones for first and second order freeness). For this, note that 
we can rewrite the crucial condition for our product of partitioned 
permutations, 

2|V| - |7r| +2|W| - \a\ =2|VV>V| - \7ra\, 

in the form 

(|V| - IttI) + - \a\) + (|V| + IWI - IVVWI) = (|VV>V| - \7ra\). 

Since all terms in brackets are non-ncgativc integers this formula can be 
used to obtain explicit solutions to our factorization problem for small 
values of the right hand side. Essentially, this tells us that factorizations 
of a disc permutation can only be of the form disc x disc; and factor- 
izations of a tunnel permutation can only be of the form disc x disc, 
disc X tunnel, and tunnel x disc. Of course, one can generalize the 
following arguments to higher order type permutations, however, the 
number of possibilities grows quite quickly. 

Proposition 5.11. 

(1) The solutions to the equation 

(l„,7n)=(0,7n) = (V,7r)-(W,(j) 
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are exactly of the form 

(1„,7„) = (0, Tt) ■ (0, TT^Sn), 

for some vr G NC{n). 
(2) The solutions to the equation 

(lm+n,7m,n) = (V, Tt) ■ (W, Cr) 

are exactly of the following three forms: 
(a) 

(lm+n,7m,n) = (0, Tt) ■ (0, TT Sm,n), 

where vr G SNc{.''^^n); 

(b) 

(lm+„,7m,n) = (0, Tt) ■ (W, TT Sm.n), 

w/iere vr G NC{rn) x NC{n) and |W| = |7r"^7m„| + 1; 

(c) 

(lm+n,7m,n) = (V, Tt) ■ (0, TT Sm,n), 

where vr G NC{m) x NC{n) and \V\ = \tt\ + 1. 

Proof. (1) The correspondence between non-crossing partitions and 
permutations was studied in detail by Biane |Bia97j . In this case we 
have 

(|V| - |7r|) + - \a\) + (|V| + IWI - |VV W|) = |1„| - \in\ = 0. 

Since all three terms in brackets are greater or equal to zero, all of 
them must vanish, i.e., 

|V| = |7r|, thus V = 0^ 

\W\ = \cr\, thus W = 0^ 

and 

|7r| + \a\ = \V\ + |W| = I VV W| = |7| = n - 1. 
(2) Now we have 

(|V|-|7r|) + (|W|-|a|) + (|V| + |W|-|VVW|) = (|VVW|-|7ra|) = 1, 

which means that two of the terms on the left-hand side must be equal 
to 0, and the other term must be equal to 1. Thus we have the following 
three possibilities. 

(a) 

|V| =|7r|, thus V = 0^, 

|W| =|(t|, thus W = 0, 

and 

|7r| + |a| = |V| + |W| = I VV W| + 1 = m + n. 
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Note that 

7rVa = VVW = 

and thus vr connects the two cycles of 7m,n- This means that vr 
is a non-crossing (m, n)-permutation. 

(b) 

|V| = |7r|, thus V = 0^, 
|W| = |a| + 1, 

and 

|V| + IWI = |V VW| = m + n- 1. 
This imphes 

|7r| + |7m,n71"~"'"| = m + n — 2, 

which means that tt must be a disconnected non-crossing (m, n)- 
annular permutation, i.e., 

TT = TTi X 7r2 with TTi G NC{m), TT2 G NC{n). 

(c) 

|V| = |7r| + l, 
|W| = |a| + 1, , thus W = 0^ 

and 

|V| + |W| = |VV W| = m + n- 1. 
This imphes 

|7r| + |7m,n71"^"'"| = TTI + Tl — 2, 

which means that vr must be a disconnected non-crossing (m, n)- 
annular permutation, i.e., 

TT = TTi X with TTi G NC{m),TT2 G NC{n). 

□ 

Example 5.12. We can now use the previous description of factoriza- 
tions of disc and tunnel permutations to write down explicit first and 
second order formulas for our convolution of multiplicative functions. 
1) In the first order case we have 

(38) (/*(;)(l„,7„) = (/*(;)(0,7n)= /(0,7r)(7(0,7r-Sn). 

■K£NC{n) 

This equation is exactly the formula for the convolution of multiplica- 
tive functions on non-crossing partitions, which is the cornerstone of 
the combinatorial description of first order freeness |NSp97j . (Note 
that vr~^7„ is in this case the Kreweras complement of tt.) 
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2) In the second order case we have 

,nj 

+ 5] (^(0' 7™,n7r-^)(7(V, tt) + /(V, n)g{0, n-'jm,n)) ■ 



TeNC{m)xNC{n) 
|V| = K|+1 



We should expect that this formula is the combinatorial key for the 
understanding of second order freeness. Ho wever, in this form it does 
not match exactly the formulas appearing in fM SS04j . Let us, however, 
for a multiplicative function / put, for vr G NC{n), 

(39) /i(7r) :=/(!„, vr) (vr G iVC(n)) 

and, for tti G NC{m) and 7r2 G NC{n), 

V>7ri X7r2. 

|v|=k|+i, 

VV(7riX7r2) = lm + n 

Note that in the definition of /2 the sum is running over all V which 
connect exactly one cycle of vti with one cycle of 112 ■ 
Then, with h = f * g, we have 

hilmAn) = Y /l(^)^l(^"^7m,n) 
7rG5]vc('n,n) 

+ (/2(7ri,7r2)^i(7rf^ X 7r2"Sm.n) 

TTi ,TT2 eNC(m) X NC{n) 

+ /i(7ri X 7r2)^2(7rfSm,vr2"Sn)). 
In this f orm we recover exactly the structure of the formula (10) from 



MSS04j . which describes second order freeness. The descriptions in 
terms of / and in terms of /2 are equivalent. Whereas / is multiplica- 
tive, /2 satisfies a kind of cocycle property. From our present perspec- 
tive the description of second (and higher) order freeness in terms of 
multiplicative functions seems more natural. In any case, we see that 
our convolution of multiplicative functions on partitioned permutations 
is a generalization of the structure underlying first and second order 
freeness. 

5.4. Zeta and Mobius function. In the definition of our convolution 
we are running over factorizations of (W, 7) into products (V, 7r)-(W, a). 
In the first order case the second factor is determined if the first factor 
is given. In the general case, however, we do not have such a uniqueness 
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of the decomposition; if we fix (V, vr) there might be different choices 
for (W, cr). For example, this situation was considered in Proposition 
15.111 in the case (2b). However, in the case when (W, o") is a disc 
permutation, it must be of the form (0,^-17, 7r~^7) and is thus uniquely 
determined. Note that factorizations of such a special form appear in 
our formula and thus deserve special attention. 

Notation 5.13. Let (^,7) G VS be a fixed partitioned permutation. 
We say that (V, vr) G VS is (U,'y) -non- crossing if 



The set of {U, 7)-non-crossing partitioned permutations will be denoted 
by VSiycip^,!), see Remark [5.81 

To justify this notation we point out that (1„, 7„)-non-crossing parti- 
tioned permutations can be identified with non-crossing permutations; 
to be precise 



Furthermore, 

U {(V, VTi X 712) I VTi G NC{m), 712 G NC{n), V > vr, |V| = |7r| + 1}. 

We can now also use a special multiplicative function, which we will call 
Zeta-function (, to single out such factorizations. It will be useful to 
be able to invert formula (jH^j) , which means we need also the inverse of 
( under our convolution. This inverse, called the Mobius-function /i, is 
a key object in the theory and contains a lot of important information. 

Notation 5.14. 

(1) The Zeta-function ( is the multiplicative function on VS which 
is determined by 



(2) The Mohius function fi is the inverse of ( under convolution, 
i.e., it is determined by 



(V,7r)-(0,7r-S) = (W,7). 



VSNciln,ln) = {(0.,7r) | TT G NC{n)}. 




1 if (1„, vr) is a disc permutation, i.e., if 1„ = 0. 
otherwise. 



6 = j2* (. 



Note that in general 
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It is also quite easy to see that the Mobius function exists and is 
uniquely determined as the inverse of the Zeta-function — the de- 
termining equations can be solved recursively. Indeed letting //„ = 
)u(ln,7n) and = /^(Im+n, 7m,n) we have 

= A«i,i + A«2 

= /il,2 + 2/ii/Xi,i + 2/i3 + 2/ii/X2 
= Il2,2 + 4//l//2,l + 4//i//l,l + 4//4 + 8//1//3 + 2/^2 + 4/Xi//2 
= //1,3 + 3//i//2,l + 3//2A^l,l + 3//4 + 6//i^3 + 3/^2 + 3//i//2 
= /X2,3 + 2/^1/^1,3 + 3/ii;U2,2 + 3/i2/Ul,2 + 9/ii/ii,2 + QfJ'lfJ'2fJ'l,l + 6/^5/^1,1 

+ 6//5 + 18yUi/i4 + 12//2Ai3 + I8//1//3 + I2//1//2 + 6/X^//2 

= jU3,3 + 6jUl//2,3 + 6a*2A*1,3 + + ^14/^2,2 + 18a*i//2jWi,2 + 18jU?jUl,2 

+ 9/^2/^1,1 + 18/ii|U2/ii,i + 9/ii/ii,i + 9/^6 + 36/UiyU5 + 27/X2/U4 + 54/^,^/^4 

+ 9/^3 + 72/xi//2A*3 + 36/X1//3 + 12/^2 + 36/^1/^2 + 9/^1/^2 

This shows how, knowing the first order Mobius function the second 
order Mobius function ^m,n can be calculated recursively. 
One should observe that with these notations we have 

In the following we will use the notation 

p-times 

It is clear, by definition, that C*^ counts factorizations into the product 
of p disc permutations, thus we have the following result. 

Proposition 5.15. For (W, 7) G VS and p > 1 we have 

C*^(W, 7) := #{(7ri, . . . , Tip) I (W, 7) = (0, TTi) ■ ■ ■ (0, TT,)}. 

Of special interest for us is the case p = 2. 

Proposition 5.16. We have for all r > 1 and n{l) , . . . , n{r) e N, 
n :— n{l) + • • • + n{r) that 

(C * C)(ln, 7n(l),...,n(r)) = #5jvc(ri(l), . . . , n{r)). 
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Proof. As noted above, (C*C)(ln5 7n(i), -- ,n(r)) counts the number of fac- 
torizations of (In, 7n(i),...,n(r)) i^to a product of two disc permutations, 
i.e., the number of factorizations of the form 

(ln,7n(l),...,n(r)) = (0, Tt) ■ (0, 7r"^7n(i),...,n(r) ) , 

with 

|7r| + |7r~"^7| = = n — r 

and vr V 7 = In. But this describes exactly connected . . . , n(r))- 

annular permutations vr G SNc{n{l), . . . , n(r)). □ 

Notation 5.17. We put 

c„(i),...,n(r) := #S'Arc(ri(l), . . . , n(r)). 

Note in particular that c„ counts the number of non-crossing parti- 
tions of n elements and thus is the Catalan number 

1 /2n 



n + l\n 

and that Cm,n counts the number of non-crossing (m, n)-annular per- 
mutations, and thus |MN04j 

_ 2mn (2m - 1\ (2n - 1 

Cm,n = ; I I I 

m + n\ m J \ n 

More generally, an explicit formula for the number of factorizations 
into p factors was derived by Bousquet-Melou and Schaeffer [BMSQOj , 
namely one has (with n := n{l) + ■ ■ ■ + n{r)) 



c--(i„a.,,,...,„,.,)-p j!;^;>!;;i;j, n[»(o( 

and thus in particular 



^pn{i) — 1 
n(i) 



Cn(l),...,n(r)=2— — — — II 



n-r + 2)\ 



n\i] 



2n{%) - 1 
nil) 



For our purposes, however, the following recursive formula for the 
number of factorizations is more interesting. 

In the next theorem we will show how to reduce the problem of 
counting the number of disc factorizations on [n] to counting the fac- 
torizations on [n — 1] . This will enable of to obtain a recursive formula 
for c„ 



-^ni,...,nr • 
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Notation 5.18. Let (W,7) be a partitioned permutation of [n] witli 
7(1) 7^ 1. Let 7jfc be the restriction of (1, A;)7(l, 7~^(/c)) to the invariant 
subset [2, n] := {2, 3, 4, ... , n}. Then 

I7I if 1 and k are in different cycles of 7, 
I7I — 1 if A; = 1 or 7(1) 
I7I — 2 if 1 and k are in the same cycle of 7, 
but A; 7^ 1 and A; 7^7(1), 

Let U — W| [2,n] be the restriction of U to [2, n] , i. e. if the blocks of U are 
Ui,. . . ,Ur and 1 G Ui, then the blocks of U are C/i, C/2, • • • , f/r where 
Ui = Uin[2,n]. In the theorem below we sum over a set of partitions 
Vk of [2, n] described as follows. 

For A; = 1, 7(1) or k not in the 7-orbit of 1, Vk = {W} i-e. Vk 
consists of the single partition U. 

For k in the 7-orbit of 1 but A; 7^ 1, 7(1), Vk = {U \ ik < |W| = 
\V(\ — 2, and U —U\/{k, 7~^(A;))}. In words this means Ui is split into 
two blocks: 

o the first containing the cycle of 7^ containing 7"^ (A;) and some (pos- 
sibly none) of the other cycles of 7 contained in Ui 

o the second containing the cycle of 7^ containing k and the remaining 
(possibly none) cycles of 7 contained in [/i. 

More explicitly, in the case k is in the 7-orbit of 1 but A; 7^ 1, 
7(1), let us write 7 as as a product of cycles di---ds where di = 
(1,7(1), . . . ,7*(1)) is the cycle that contains 1. Let d[ = (7(1), 7^(1), 
. . . , 7~^(A;)) and < = (A;, . . . , 7*(1)). Then % = d[d'ld2 ■■■ds. Vk con- 
sists of all partitions U of [2,n] such that U = {U[,Ui ,1/2, ■ ■ ■ ,Ur} 
where U[ U ?7(' = Ui,U[n U'l = 0, U[ contains d[, U'l contains <, and 
each cycle of 7 that was in Ui is now in either t/{ or t/f, i.e. % < U 
and |W| = |W| - 2. 

Theorem 5.19. 

(41) c\un)-j2T.C\ii,fk) 

Proof. We must show that for each factorization (0, tt) ■ (0, a) of (W, 7) 

there are k := 7r(l), U G Vk, and permutations of [2,n], n and a such 
that (0,71") ■ (0,(7) = (V(,^k)- Conversely we must show that given k, 
U E Vk and a factorization (0, tt) ■ (0,0") of {U,%) there are n and a 
such that (0, tt) • (0, a) = {U, 7) and 7r(l) = A;. Moreover we must show 
that these two maps are inverses of each other. The relation between 
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TT, cr and TT, a is given by tt = (1, A;)7r|[2,n], a = cr(l,7 ^{k))\[2,n]- So on 
the level of permutations we have a bijection. The main work of the 
proof is to show that starting with tt and a we have U :— n V a E Vk 
and 2\U\ — \^k\ = l^r] + \a\; and then conversely starting with hi E Vk 
and a factorization (0, tt) • (0, o") of {IA,%) then 2|W| — |7| = |7r| + \(t\ 
and Tiy a = U. 

Note that we have for all k 

fkl-l A; 
' ' [|7r| k^l 

1^, fkl-1 A: ^7(1) 
\\a\ fc = 7(l) 

It is necessary to break the proof into four cases: k is not in the 
7-orbit of 1; k is in the 7-orbit of 1 but k ^ 1,7(1); k = 1; and 

A; = 7(1)- 

Suppose we have a factorization 

(W,7) = (0,7r)-(0,(j), 

i.e., 7 = TTcr, U — TT V a, and 

2\U\ - I7I = |7r| + \a\ 

with A; := 7r(l) not in the 7-orbit of 1. Then|7fe| = I7I and Vk contains 
only the partition of [2,n] which results from U by removing 1, i.e. 
lA = IA. Then we have \IA\ = — 1. Hence |7r| + l^l — \tt\ + |cr| — 2 = 
2\U\-\^\~2=\U\^\-f\ = \U\-\%\. 

Also 0^|[2,„] = O^i- and 0^|[2,n] < 0;^^. Thus W = (tt V 7)|[2,n] < tt V 7^. 
On the other hand the difference between 0^|[2,„] and 0^^ is that the 
blocks containing 1 and k have been joined. However these points were 
already connected by tt. Thus tt V 7^ <U, and so U — it \/ a, and thus 

(W,7) = (0,7r)-(0,a). 

Conversely, given a factorization (0, tt) • (O,^) of {U,^k), let tt = 
(l,A;)7r and a = a{l,'-f~^{k)). Then n W a = U because 1 has been 
connected to the block of U containing k. Also 7^(7r) = #(71") and 
#((t) = thus |7r| = |7r| — 1 and |cr| = |a| — 1, and so |7r| + \a\ = 

2\U\ — |7|. This estabhshes the bijection when k is not in the 7-orbit 
of 1. 

Let us now consider the case that 1 and k are in the same cycle of 
7, but k 7^ 1,7(1). Again suppose that (0,7r) • (0, a) is a factorization 
of (W, 7) with 7r(l) = k. In this case we have that |7fc| = I7I — 2 and so 
by the triangle inequality. Lemma 4.7 
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2|7r V a| - I7I + 2 = 2|7r V a| - I7I 
= \{n y a, na)\ 
< |(0,7r)| + |(0,a)| 
= |7r| + l&l 
= |7r| + \a\-2 
-2|W|-|7|-2, 

and thus 

Itt V(t| < \U\ -2. 

On the other hand, let us compare 

V (J = TT V 7 with U = n V 'j. 

Note that all our changes of the permutations affected only what hap- 
pens on the first cycle of 7. Since the transition from 7 to 7 consists in 
removing the point 1 and splitting the first cycle of 7 into two cycles, 
we can lose at most one block by going over from tt V 7 to tt V 7. Thus 

Itt V o-| = (n - 1) - #(71 V 0-) > (n - 1) - (#W + 1) = |W| - 2, 

so that we necessarily have the equality 

Itt V<t| = \U\ - 2. 

Thus U := Ti\/ a E Vk and 2|7rV(j| — |7fc| = |7r| + \a\. Hence (0, tt) • (0, a) 
is a factorization of {U,jk)- 

Conversely let us suppose that k is in the 7-orbit of 1 but A; 7^ 1 or 
7(1) and U G Vk and (0, tt) • (0, a) is a factorization of {U, 7^). We must 
show that TT V (7 = W and that |7r| + \(x\ = 2\IA\ — 1 and k are in the 
same orbit of tt and 1 and ^i^^{k) are in the same orbit of a. So the 
blocks of lA containing d'^ and ri" arc joined in tt V o". Thus n \/ a = U. 
Also \U\ = \U\ - 2, so |7r| + \a\ = |7r| + |<t| + 2 = 2\U\ - \^k\ + 2 = 
1^1 ~ l7fc| ~2 = 2\U\ — I7I. Thus (0, tt) • (0, a) is a factorization of (U, 7). 
This establishes the bijection in the case k is in the 7-orbit of 1 but 
A; 7^ 1 or 7(1). 

Next suppose that k — 1 and (0, tt) • (0, a) is a factorization of (W, 7) 
with 7r(l) = 1. Then |7r| + |<t| = |7r| + |cr| - 1 = 2\U\ - I7I - 1 = 
2\U\ — I7I -|- 1 = 2\U\ — |7fc|. Let Ui be the block ofU containing 1 and 
Ui = Ui (1 [2,n]. We must show that f/i is a block of tt V 7^. Since 
Try ■J = U we know that if di and dj are cycles of 7 contained inUi then 
TT must connect them. Since Tr\jj^ — Tt\jj^ we see that tt connects the 
corresponding cycles of % (which are unchanged except for the cycle 
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containing 1). Similarly if fi and /2 are cycles of tt contained in Ui 
and neither is a singleton then they are connected by 7 and thus by 
7fc. Thus (0,71") ■ (0,(3") is a factorization of (W,7fc). 

Conversely suppose that k = 1, U G Vk, and (0, tt) ■ (0,(3") is a 
factorization of (W, 7^). We must show that 7r(l) = 1 and (0, tt) ■ (0, a) 
is a factorization of (^,7). Since tt V 7^ = W and 7 connects 1 to 
7(1) G f/i, we have that vr V 7 = W. Also |7r| + |(t| = |7r| + |(3"| + 1 = 
2|W| - |7fc| + 1 = 2\U\ - |7i| - 1 = 2\U\ - |7|. Thus (0,7r) ■ (0, a) is a 
factorization of (W, 7). This completes the case when k = 1. The proof 
in the case k = 7(1) is exactly the same except that the roles of tt and 
a are reversed. □ 

Let us take a closer look at the meaning of Theorem 15 . 1 91 for the case 
(W,7) = (In, 7n(i),...,n(r))- To rcducc the depth of subscripts we shall 
write c(ni, . . . , n^) for c„(i),...,n(T-)- 

Proposition 5.20. We have for a// r, ni, . . . , ri^ G N the recursion 

r 

(42) c(ni, . . . ,nr) = ^ni ■ c(ni + - 1, ns, . . . , n;_i, n^+i, . . . , n^) 

1=2 

+ Y1 c{k-l,ni^,...,nijc{ni-k,nj^,...,nj,) 

k=l A={ii,...,is} 
B={ji,...,jt} 



where the sum is over all pairs of subsets A,B G [2, r] such that ACiB = 
and AU B = [2, r] including the possibility that either A or B could 
be empty. We have for all m,n > 1 

l<fc<n 

and 

(43) Cj7i,n ^ ^ {c}i—\Ciji,n—k Ciji,k—lCn—k) ^Cm+n— 1, 

l<k<n 

where we use the convention that cq = 1 but c(ni, . . . ,nr) = if r > 1 
and for some i, rii = 0. 



Proof. Let n = riiH hn,.. By Proposition|5J^c(ni , . . . ,nr ) = C (1 

7ni,...,nr)- So we must give the correspondence between the terms on 
the right hand side of (j4H) and the right hand side of (j42|) . In this case 
U = In and U = l^-i (in the notation of I5.18|) . Thus Vk = {l„_i}. 
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Also for ni H — ■ + n/_i < k < ni + ■ — \-ni, (*^{ln-i, 7a:) = c{ni +ni — 
1, n2, . . . , ni^i, m+i, ...,nr). Thus 

n n 

(44) E Y.c\u.%)= E c\in-,.ik) 

fc=ni+l^g-p^ fc=ni+l 
r niH hri; 

= E E c(ni + - l,n2, . . . ,ni_i,r2;+i, . . . ,72^) 

/=2 fc=niH hrii— 1+1 

= ^ni- c{ni +ni-l,n2,..., ni_i, m+i, ...,nr) 

1=2 

For k < ni, jk = d'^d'{d2 ■ ■ ■ dr, with d[ a cycle of length k — 1 and c?'/ a 
cycle of length rii — k. Vk is the set of all partitions of the cycles of 7^ 
into two blocks such that d'l and d'l are in different blocks. Hence 

EC*'(^,7fc)= E c(ni-/c,ni,,...,nijc(/i;-l,njj,...,njj 

ueVk A={ii,...,is} 

B={ji,...,jt} 

where the sum is over all pairs of subsets A,B G [2, r] such that ACiB = 
and AU B = [2,r] including the possibility that either A ot B could 
be empty. Thus 

rii 

(45) EEC*' (^,7.) 

ni 



E E c{ni-k,ni,,...,nijc{k-l 



1 "^ji 1 • • • 1 "^jt ; 



k=l A={ii,...,is} 
B={ji,...,jt} 

Assembling equations and (ji^ gives the result. □ 

In |QZ84j ■ O'Brien and Zuber used a similar formula of this kind in 
order to compute the asymptotics of, so called, external field matrix 
integral. See also |BMS00j and Theorem 15.221 

Clearly, our notions around the convolution of functions on VS are 
analogous to (and motivated by) the convolution of functions on posets. 
Even though we are not able to put the above theory into the frame- 
work of posets, it seems that this analogy goes quite far. The following 
description of the Mobius functions is an instance of this — its poset 
analogue is due to Hall (see |Rot64j ). It is essentially the simple obser- 
vation that one can expand the Mobius function in terms of a geometric 
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series as 

oo 

Proposition 5.21. We have for any (^,7) G VS that 

00 

MW,7) = 5(W,7) + E E (-l)^ 

k=l (W,'y) = {0,7ri)...(0,7rfc) 

Proof. As noted above this is just the geometric series for 

(s+{c-5)r-\ 

(Note that we are working for this in the algebra of functions on VS 
with the pointwise sum and the convolution as sum and product — we 
are not bothering about multiplicativity.) The only thing to check is 

that the sum is finite, and this is the case because the number of factors 
k is bounded by \[U,'y)\, since |(0,7r)| > 1 for any tt 7^ e. □ 

This description of the Mobius function allows us now to derive a 
recursive formula for jj,. 

Theorem 5.22. Consider (^,7) G VS such that 7(1) ^ 1. Then we 

have 

(46) M^^,7) = (-l) Yl MV,7r), 

(0,(l,fc))-(V,^)=(M,7) 

where the sum runs over all decompositions of {U, 7) into a product of 
a disc transposition (0, (1, k)) (with k > 2) and a (V, tt) G VS. 

The proof of this theorem will rely on the following lemma. 

Lemma 5.23. Let {Uri) e VS such that 7(1) ^ 1. For p G N, we 
denote by Sp the set consisting of all tuples (tti, . . . , tt^) of permutations 
such that TTi ^ e for all i — 1, . . . ,p and 

(0,7ri)---(0,7rp) = (W,7). 

We consider now the two sums 

00 

(47) S,:^Y1 E (-1)' 

P=l (7ri,...,7rp)e5p 

and 

00 

(48) S^-.^Yl H (-1)' 

p=l (iri,---,Tp)eSp 

TTI = (1, k) for k^l 
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where the second sum S2 is over all tuples (tti, . . . , VTp) as for the first 
sum Si, but now with the additional property that tci is a transposition 
interchanging the element 1 with some other element. 
Then the two sums and PHj) are equal, 

Si = 5*2. 

Proof. Let tt = (tti, . . . , VTp) G Sp. Let 1 < q < p denote the smallest 
index for which 1 is not a fixed point of vr^ ; note that such a q necessarily 
exists since 7(1) 7^ 1. We shall group all factorizations into three 
classes: la), lb) and 2). Class 1) consists of factorizations for which 
7Tq is a transposition interchanging 1 with some other element. The 
subclass la) consists of factorizations for which q = 1 and subclass lb) 
of those for which q > 2. Class 2) consists of all other factorizations. 
Let n = (vTi, . . . , TTp) be a factorization from the class lb). We define 

W = . . . , 7rp_ J = (tTi, . . . , 7rg_2, TTg^iTTg, TT^+i, . . . , TTp). 

In the following we shall prove that f : U \—>- U' is a bijection between 
factorizations of class 16) and factorizations of class 2). 

Firstly, we prove that U' G Sp and is of class 2). Clearly, 7r^_^ = 
Tiq^iTTq IS a pcrmutatlou which does not fix 1, it is not a transposition 
interchanging 1 with some other element, and we have 

(0,7r,_i)-(0,7r,) = (0,7r;_i). 

In order to show that / is a bijection we shall describe its inverse. If 
n' = {tt[, . . . , 'n'p_i) E Sp and is of class 2), we define l<g<p— Ito 
be the smallest number for which nq_i does not fix 1. There is a unique 
decomposition 'n''q_i = Tiq-ii^q such that 1 is a fixed point of TTq-i and 
TTg is a transposition interchanging 1 with some other element. Thus 
|7rg_i| + iTTgl = |vrg_]^|. The assumption that the factorization 11' is of 
class 2) implies that Hq-i 7^ e. For 1 < i < g — 2 we set tt; = tt^' and for 
g + 1 < i < p we set TTj = 'n:[__i. In this way we defined 11 = (tti, . . . , vr^). 
Now it is easy to check that : 11' i— > 11 is a left and right inverse of /. 

Since the factorization 11 and the corresponding 11' contribute to (j47|) 
with the opposite signs, the contribution of all factorizations of class 
lb) cancels with the contribution of factorizations of class 2). □ 

Proof of \5.2^ In the proof we will consider all factorizations (0, tti) ■ 
(0, 712) ■ ■ ■ (0, Hp) = (U, 7) with the requirement that vTj 7^ e for all 
i.e. (vTi, . . . , TTp) G Sp, as in the proof of Lemma 15.231 Sometimes we 
will require in addition that tti = (1,A;) with k ^ 1. To simplify the 
notation we will not explicitly state every time that tTj 7^ e. Since 
7(1) 7^ 1 we have 6(U,'j) = 0. When 7 is a transposition the right 
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hand side of equation (pUj) is —1; so we can assume that 7 is not a 
transposition. So by Proposition 15.211 we have 



M^^,7) = E E i-ir^Y: E (-1) 



00 



P=l (0,7ri)---{0,7rp) P=l (0,{l,fe))-(0,7rp) 

={W,7) =(2^,7) 

00 

= E E E (-1)' 

P=2 {0,(l,fe)),(V,7r) (0,7r2)-(0,7rp) 
(0,{l,fc))-(V,7r) = (W,7) ={V,7r) 
00 

= - E EE (-ir' 

(0,{l,fc)),(V,7r) p=2 {0,7r2)---(0,7rp) 
{0,{l,fc))-(V,7r) = (W,7) =(V,7r) 
00 

= - E EE i-ir" 

(0,{l,fc)),(V,7r) p=2 (0,7r2)---(0,7rp) 
{0,{l,fc))-(V,7r) = (W,7) =(V,7r) 



(0,{l,fe)),(V,7r) 
{0,{l,fc)).(V,7r) = (W,7) 



□ 



One observes that the recursion formulas for the Mobius function 
and for (^*^ look very similar. However, there are some significant dif- 
ferences. The recursion for (*'^ effectively expresses C*^ for ^ points in 
terms of C*^ for n — 1 points. The recursion for the Mobius function 
does not reduce the number of points. Nevertheless, at least for first 
and second order one can match the two recursions and connect the 
values of the Mobius function with the values of the function (*'^ (i.e., 
with the number of non-crossing partitions and non-crossing annular 
permutations). In order to see this let us first specify the meaning of 
Theorem 15.221 for first and second order. In first order we get 

/^(ln,7n)=- E /"(lfc'7fe)/^(ln-fe,7n-fc), 
l<fc<n-l 

which shows that (— 7^+1) and C*^(ln,7n) satisfy the same 
recursion (namely the one for the Catalan numbers). This is, of course, 
just the well-known fact |Kre72t Spe94| that the Mobius function on 
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non-crossing partitions is given by the signed and shifted Catalan num- 
bers. In second order our recursion reads 

k)l^{ln-k,'yn-k) + k)l^{h,1k)), 

l<fc<n-l 

which we recognize — by taking into account the shifted relation be- 
tween ij, and C*^ on the first level — as the recursion for (— 
1m,n)- Let US collect these explicit results about the Mobius function 
in the following theorem. 

Theorem 5.24. We have for m,n eN that 

7„) = (-1)"-^ • #7VC(n - 1) = (-1)"-^ • c„_i 

and 

For higher orders we were not able to match the values of /i with 
those of C*^. 

6. R-TRANSFORM FORMULAS 

Let us consider the situation that two multiplicative functions / and 
h on VS are related hy h — / * C- We want to understand what 
this means for the relations between the numbers /«„ := /(1„,7„) and 
Km,„ := /(lm+n,7m,n) ou ouc sidc and the numbers a„ := /i(l„,7„) 
and am,n '■= h{lm+n,'lm,n) on the other side. In particular, we want to 
express this in terms of the generating power series of these numbers, 

C{x) := 1 + XI '^nx'', C{x, y) := ^ Km,nX"Y 

n>l m,n>l 

and 

n>l m,n>l 

(Note that the above summation corresponds to putting formally 

/(lo,7o) := 1 and /(lo,7o,o) := 

for a multiplicative /. Our notation is motivated by the fact that the 
most important realization of the relation h = f*C, will be the situation 
where the a's are the correlation moments and the k's the correspond- 
ing cumulants, thus M is a moment series and C is a cumulant series.) 
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On the first order level we have 

TT(^NC{n) 

which is the usual moment-cumulant formula of free probability theory, 
and it is well-known |Spe94| that this is equivalent to 

C{xM{x)) = M{x). 

Our main goal now is to derive the analogue of this for the second 
order level. There we have 

'reS'jvc(''n,n) ttj X7r2eiVC(m)xJVC(n) 

V| — |7r]_ X7r2l+1 

It turns out that the second term, the sum over disconnected par- 
titions, is quite easy to deal with. The first term, the sum over con- 
nected annular permutations, looks much more involved, however, one 
can handle this also if one realizes that one can reduce this first term to 
the second one. Namely, one can sum over all connected annular permu- 
tations by first bundling all through-cycles into one through-cycle and 
secondly decomposing this through-cycle into sub-cycles all of which 
are through-cycles. In this way one can reduce the problem of dealing 
with all annular non-crossing permutations to the problem of consid- 
ering permutations with exactly one through-cycle and the problem 
of considering permutations where all cycles are through- cycles. The 
first problem corresponds exactly to the above sum over disconnected 
partitions. So we can write 



Y /(0-^)= Y /(V,vriX7r2), 

■K&Sl^lc{rn,n) ttj xirj eiVC(m) x ]VC(n) 

V| — |7r]_ X7r2|+1 

where / is now the multiplicative function corresponding to 

f iX'fii In) ^ri; /'(Im+ri; '~^m,n) f^m,n 

with 
and 

7r€5^'L(m,n) 
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Thus we can combine this to get finally 



J2 (/(V,7riX7r2) + /(V,7riX7r2)) 

Tl x-n-2 eiVC{m) X NC(n) 
|V| = kiX7r2l + l 



TTi X TTQeNCim) X NC(n) 
|V| = |7riX7r2|+l 

where g is the multiplicative function corresponding to 

fl'(ln)O'n) Oiji, g(^ljTi-^nT 'ym,n) Q^m,n 

with 
and 

So we have to translate the relation between K.rn,n and / and the relation 
between a^.n a-nd into relations between the corresponding formal 
power series. 

Proposition 6.1. Let f be a multiplicative function on VS with 
/(!„, 7„) =: Kn and C{x) := 1 + ^ KnX^. 

n>l 

Put 

where S^Q{m,n) denotes the permutations in SjsfciiTijn) for which all 
cycles are through- cycles. Consider the corresponding generating power 
series 

C{x,y) := J2 i^m,nx^y''- 

m,n>l 

Then we have 

Cix v) - -x^^loef^^M^^l 
C{x,y)- xy^^^^log[ ), 

or equivalently 

f {C{x)-xC'{x)){C{y)-yC'{y)) 1 ^ 

^"''^"""''^ {xC{y)-yC{x)Y (^^^ 
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Proof. Note that we can parametrize an element tt e 5'^[;(m, n) in a 
bijective way by specifying the number of cycles, the number of ele- 
ments on each circle for all cycles, the position of a fixed element (let's 
say 1) in its cycle and the first element on the other circle of this cycle. 
Let us denote the number of cycles by r, the number of elements of 
the cycles on the first circle by ii, . . . ,ir and the number of elements 
of those cycles on the other circle by ji, . . • , j, • Thus the l-th cycle 
contains ii +ji elements and makes the contribution in the calcu- 

lation of Km,n- We normalize things so that the first cycle contains the 
element 1. Fixing ii, . . . , and ji, . . . ,jr we thus have ii possibilities 
for where 1 sits in the first cycle and n possibilities for the first element 
of this cycle on the other circle. This means we have 



nH Vir=m JiH \-jr=n 



and thus 

r>l n,...,ir>l jl,...,jr>l 
r>l il,...,ir>l jl, — ,jr>l 



d_ 

r>i jl,jl>l i2,j2>l V,jr>l 



Yy^{i Y Hf^n+ji^^'y'') ■ i Y ^^^2+n^''y'')■■■[ Y ^^^r+oyy'')) 



dy 

Let us now use the notation 



c{x,y) ■.= Y i^i+J^y- 



Then we can continue with 

d d 



r>i y 
= ^y^§^{Ylc{x,yy) 

^ r>l 
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The assertions follow now by noting that 

xC{y) - yC{x) 



x-y 

and by working out the partial derivatives. □ 
Proposition 6.2. Let g be a multiplicative function on VS. Put 

and denote its generating power series of second order by 

H{x,y) := J2 ^m,nx^y''- 

m,n>l 

Put 

ttn := (5'*C)(ln,7n) 

and 



{V,-K-l X772) 

|V| = |7riX7r2l + l 



and denote the corresponding generating functions by 



M{x) := 1 + ^ cKnx" and M{x, y) := ^ a^^nx""?/". 

n>l m,n>l 

Then we have the relation 

M{x,y) = H{xM{x),yM{y)).{l + x^) ■ + 

Proof. Let us do the summation in the definition of „ in the way 

that wc first fix the two cycles Vi G tti and V2 G which are con- 
nected by V and sum over all possibilities for fixed Vi,V2. If Vi has 
k elements and V2 has / elements then this contributes the factor ak,i. 
Furthermore, 7ri\Vi decomposes into k independent non-crossing par- 
titions and the summations over them (for fixed Vi) gives the ccj for 
the intervals between consecutive elements from Vi. (Of course, we are 
counting here modulo m.) For the final summation over Vi we have to 
notice that there are two different possibilities: either a fixed number 
(let's say 1) is an element of Vi - in which case we can specify the situ- 
ation by prescribing the number k of elements of Vi and the differences 
ii, . . . ,ik between consecutive elements in - or 1 is not an element of 
Vi, — in which case we need an extra factor ii, because we have now ii 
different possibilities how 1 can lie between two consecutive elements 
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of Vi. Since we have the same situation for V2 we can thus write a^, 
in the form 



a. 



k,l>l <i,--,ifc>o h:--,n>o 

fc+ilH +ii^=m i+JiH i-j(=n 



Translating this into generating power series gives the assertion. □ 
The combination of the previous two propositions, with 
H{x,y) = C{x,y) + C{x,y), 
gives now our main result. 

Theorem 6.3. Let f and h be multiplicative functions on VS which 
are related by 

Denote 

and 

and define the corresponding generating power series 
C{x) := 1 + ^ C{x, y) := ^ 

n>l m,n>l 

and 

M(x) := 1 + ^nx"", M(x, y) := J] 

n>l m,n>l 

Then we have as formal power series the first order relation 

(49) C{xM{x)) = M{x) 
and for the second order 

(50) M{x,y) = H{xMix),yMiy)) • ^^j^ ■ 
where 

(51, «(..,)..C(..„-..4log(£^M^). 
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or equivalently, 

±{xM{x)) i{yM{y)) 



(52) M{x,y) = C{xM{x),yM{y)) 

+ xy 



M{x) M{y) 



" 'xM{x)) ■ j-iyM{y)) 



{xM{x) -yM{y)Y {x - yf . 

Proof. The formulation ()50|) and ()5H) follows directly from a combina- 
tion of Propositions 16.11 and 16. 2[ In order to reformulate this to 
one uses the equivalence of the two formulas in Proposition 16. II and the 
fact that C{xM{x)) = M{x) yields 

1 - xC\xM{x)) = ,f^[f, 

£{xM{x)) 

□ 

If we go over from the moment generating series M to a kind of 
Cauchy transform like quantity G, then these formulas take on a par- 
ticularly nice form. 

Corollary 6.4. Consider the same situation and notations as in The- 
orem \6.3i . In terms of 

G{x) := -M(l/x), G{x, y) := — M(l/x, 7^(x, y) := — C(x, y) 

X xy xy 

the Equation \5^) can he written as 
(53) 

G{x,y) = G\x)G\y){n{G{x),G{y)) ^ ^ ^ 



{G{x)-G{y)Y} {x-yy 

TZ{x, y) is the second order i?-transform. Note that Voiculescu's first 
order i?-transform IZ is defined by the relation C{x) = 1 -|- z7l{x), and 
equation (jl^ says for this 

-l^ + 7^(G(x))=x, 
G{x) 

i.e., that G{x) and K{x) ■= ^ + lZ{x) are inverses of each other under 
composition. 

Example 6.5. Let us apply our formulas to some examples. 

1) If we put / to be the multiplicative function with ^2 = 1 and 
all other k„ and all Km,n vanishing, then h = f * ( counts the non- 
crossing pairings, i.e., in this case M{x) is the generating function of 
the number of non-crossing pairings (on one circle) and M{x, y) is the 
generating function of the number of non-crossing annular pairings (on 
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two circles). Let us calculate it by using the above theorem. 
We have 

C(a;) = l + x^ C{x,y) = 
and we know that M is the generating function of number of non- 
crossing pairings on a circle. In this case 

C{x,y) = xy, 

and thus 

H{^^ y) = -^y^ iog(i - xy) = 

which yields the result 

, £(xM(x)) • jiyMjy)) 

M [x, y) = xy ■ 5—. 

[I- xyM{x)M{y)) 

Related formulas are known in the physical literature, see, e.g, |FMP78j . 
jSZHH], |KKP95j . 

2) If we put / = C then h = ( * ( counts the non-crossing permu- 
tations, i.e., in this case M is the generating function of the number 
of non-crossing permutations (which is the same as non-crossing parti- 
tion) on one circle and M2 is the generating function of the number of 
annular non-crossing permutations (on two circles). 
We have 

C{x) = -^, C{x,y)=0. 
1 — X 



In this case 



r^f N l-x-y 



il-x)il-y) 



and thus 



d'^ xy 
H{x, y) = -xy—— log(l - xy) 



dxdy (1 — X — yY' 

which yields 



M(x, y) = xy 



A 



-XxM{x)) ■ j-(yM{y)) 



,2- 



(1 -xM(x) -yM{y)y 

3) Let us finally see whether we can extract the value of the Mobius 
function from our formula. Since we have 5 = * C, our formula with 

M(x) = l + x, M(x,?/) = 

should allow to solve for C{x,y) which is then the generating function 
for the annular Mobius function. Note that we already know M{x) in 
this case to be the generating function of the disc Mobius function. 
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If M{x, y) vanishes identically this implies that H{x, y) vanishes iden- 
tically, leading to the identity 



xy 



dxdy X — y 

iCix)-xC'ix))-{Ciy)-yC'iy)) 



{xC{y) - yC{x)y {x - yf 



Remark 6.6. Equation l52l gives the second order version of moment- 
cumulant relations. 

"1,1 = '«l,l + 1^2 
«2,1 = 1^1,2 + 2Ki/ti,i + 2/t3 + 2KiK2 
"2,2 = ^^2,2 + 4kiK2,i + Ak\ki^i + 4^4 + + 2k\ + Ak\k2 

(^1,3 = ^1,3 + 3KiK2.i + 3K2/tl,l + 3fi;4 + 6/ti/t3 + ?>k\ + ?)k\k2 
"2,3 = «^2,3 + 2KiKi,3 + 3KiK2,2 + 3K2Ki^2 + ^f^\f^l,2 + 6Ki/t2Kl,l + 6k^Ki_i 

+ 6^5 + I8K1K4 + I2K2K3 + 18k^K3 + 12kik\ + 6k'^K2 

«3,3 = «;3,3 + 6KlK2,3 + 6K2/tl,3 + 6K^Ki_3 + 9/t^K2,2 + 18Kl/t2Kl,2 + 18KiKi_2 
+ Qk\ki^i + 18/ti/t2/tl,l + Qk\ki^i + 9^6 + 36K1K5 + 27/t2/t4 + 54k^K4 

O Q o 0/1 

+ + 72K1K2K3 + 36k^K3 + 12k2 + 36/t^/t2 + 9K]^K2 

^^1,1 = «i — ^2 + ai,i 
Ki^2 = —Aoi\ + 6«ia2 20:3 — 2aiai 1 + ai 2 
K2,2 = 18a^ — 36a;^a2 + 6^2 + 16ai«3 ~ 4^4 + Aa\ai^i — Aaiai^2 + «2,2 
Ki^3 = 15a^— 30a^a2+6«2+-'-2«iCt3~3Q;4+6a^a;i^i— 3a2Cti,i— 3«iai^2+«i,3 
K2,3 = —12a\ + 180 02 — 72Q;ia2 ~ §4 a^a^i + 24a;2a;3 + 30«ia4 — 605 

— \2a\ai^i + Qaia20Li^i + 12a\ai^2 — 3a;2«i,2 — 2a;ia;i^3 — 3aia2,2 + 0^2,3 
fi;3,3 = 300a?-900a^a2 + 576a?a^-48a^ + 432a?a3-288aia2a3 + 18a3 

— 180aja4 + A5a2a4 + bAaia^ — 9aQ + 36q;^q;i,i — 36a;Ja2«i,i + Qa^ai^i 

— 36afai_2 + 18aia2«i,2 + 12Q;^ai_3 ~ 6Q;2ai,3 + 9a^a2,2 ~ 6q;iQ;2,3 + 03,3 



60 b. collins, j. a. mingo, p. sniady, and r. speicher 

7. Higher order freeness and corresponding cumulants 
7.1. Abstract framework. 

Definition 7.1. A higher-order (non-commutative) probability space, 
or briefly HOPS, {A, ip) consists of a unital algebra A and a collection 
= {fn)neN of maps {n G N) 

(fin '■ Ax ■ - X A C, 

n times 

which are linear and tracial in each of its n arguments and which are 
symmetric under exchange of its n arguments and which satisfy 

and 

02, . . . , a„) = 
for all n > 2 and all a2, ■ ■ ■ ,an & A. 

Of course, we can include the usual (first order) non-commutative 
probability space (.4, ipi) into this framework by putting all higher 
(fn equal to zero. In the same way we recover a second order non- 
commutative probability space {A,(pi,(p2) by putting = for all 
n > 3. 

Definition 7.2. 1) We denote by VS{A) the set of partitioned per- 
mutations decorated with elements from A, i.e., 

PS= \J{VS{n) xA""). 

2) For a function 

/ : VS{A) C 
(V, tt) X (ai, . . . , a„) /(V, 7r)[ai, . . . , a J 
and a function 

g:VS 

we define their convolution 

f*g:VS{A)^C 

by 

(/*^)(Z^,7)K,...,a„] := f{V,7r)[ai,...,an]- g{yV,a) 

(V,Tr),(W,a)eVS(n) 
(V,7r)-(>V,CT)=(W,7) 

for all (W, 7) e VS{n) and all oi, . . . , a„ e ^. 
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Definition 7.3. A function / : VS{A) C is called multiplicative if 
we have 

/(V,7r)[ai, . . . ,a„] = /(l^, 7r|s)[(ai, . . . ,a„)s] 
Bev 

and 

/(l„,o-"Vo-)[a^(i), . . . ,a^(„)] = /(l„,7r)[ai, . . . , a„] 
for all Oi, . . . , a„ e .A and all tt, cr e S{n). 

Note that this extension of our formalism on multiplicative functions 
on VS and their convolution from the last section is not changing the 
results from the last section. The structure of all formulas remains the 
same; one just has to insert the dummy variables at the 

right positions. Thus, in particular, S is still the unit for this extended 
convolution and f = g * ( is equivalent to g = f * fi for multiplicative 
/, g on VS{A). And again, the convolution of a multiplicative function 
on VS{A) with a multiplicative function on VS gives a multiplicative 
function on VS{A). 

It is clear that a multiplicative function / on VS{A) is uniquely 
determined by the values of 7n(i),...,n(r))[oi, . . . , a„] (where we put 
n := n(l) + • • • + n(r)) for all r G N, all . . . ,n(r) G N and all 
ai, . . . ,an & A. 

7.2. Moment and cumulant functions. Let us now apply this for- 
malism to get moment and cumulant functions for higher order prob- 
ability spaces. So let a HOPS {A, (f) be given. We will use the (pn to 
produce a multiplicative "moment" function on VS (A) , which we will 
also denote by (p. Namely, we put 

<^(ln,7n(l),...,n(r))[ai,--- 

:= (Pr{o,l ■ ■ ■ fln(l); • • • ; (ln{l)+-n{r-l)+l " " " dn) 

and extend this by multiplicativity. (Note that we need the </?„ to be 
tracial in their arguments for this extension.) 
Here is an example for our function ip. 

(^({1,3,4}{2}, (l,3)(2)(4))[ai, 02, 03,04] = (^2(0103,04) • <^i(a2) 

Definition 7.4. For a given HOPS {A, </?) we define the corresponding 
(higher order) free cumulants as a function on VS{A) by 



K — (p * fJ,, 
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or more explicitly 

K(W,7)[ai,...,a„] (^(V, 7r)[ai, . . . , a„] • /x(>V, cr), 

(V,7r),(W,<T)SP5(n) 
(V,7r)-(W,a-)=(W,7) 

for all neN, {U, 7) G VS{n), Oi, . . . , a„ e ^. 

As we noted before the definition above is equivalent to the statement 
= K*C, i-e., 

ip{U,'y)[ai,...,an\= k(V, 7r)[ai, . . . , a„] 

(V,7r)e7'5jvc(W,7) 

for all (W,7)[ai,...,a„] G 

Furthermore, as with ip, k is also a multiplicative function on VS{A). 
Thus in the same way as all ip{U,'y) are determined by the knowledge 
of all 

¥'(ln,7n(l),...,n(r))[ai, • • • , ^n] 

= ■ ■ ■ an(l); • • • ; '^n(l)+--+n(r-l)+l " ' " (^n{l)+-+n{r)) 

the free cumulants k{U, 7) are determined by the values of 

«(ln,7n(l),...,n(r))[ai, . . . 

= • '*n(l),...,n(r)(cH) • • • ) 0,n{l)] ■ ■ ■ ] Ctn(l)H |-n(r-l)+l) • • • ) Ctn(l)H |-n(r))- 

Remark 7.5. Note that whereas on the level of (p we also know (by 
definition) that wc can multiply elements along the cycles of vr (and 
thus we do not need a comma as separator for those elements along a 
cycle), this is not true for k. Thus we have, e.g., 

V^ih: (I,2)(3))[ai,a2,a3] = 932(0102; 03) = ^{h: (1), (2))[aia2; 03], 

but no clear relation exists among 

k{13, (1, 2)(3))[ai, a2, 03] = «;2,i(ai, ^2; 03) 

and 

fi;(l2, (1), (2))[aia2; 03] = Ki,i[aia2; 03]. 

Note also that since our convolution on VS coincides on the first level 
with the usual convolution of multiplicative functions on non-crossing 
partitions, the above definition of cumulants reduces on the first level 
to the usual free cumulants. 
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7.3. Higher order freeness. Equipped with the notion of cumulants 
we can now define "freeness" by the requirement of vanishing of mixed 
cumulants. 

Definition 7.6. We say that a family (A'j)jg/ of subsets of A is free 
( of all orders ) if we have the following vanishing of mixed cumulants: 
For all n > 2 and all G Xi(k) ^ k < n) such that i{p) ^ i{q) for 
some 1 < p, g < n we have 



for all TT G S{n). 

Example 7.7. Let us see that this definition includes the definition 
of Voiculescu |VDN92j for (first order) freeness and the definition of 
Mingo and Speicher |MS04j for second order freeness. 

1) On the first level this follows from the fact that our cumulants 
reduce then to the usual free cumulants and it is well-known that free- 
ness is equivalent to the vanishing of mixed cumulants. One can see it 
directly as follows: Let us consider G Xi(^k) with i{k) ^ i{k + 1) and 
'■PiiO'k) = for all = 1, . . . , n. Then we have 

V9i(ai • • -a^) = V5(l„,7„)[ai, . . . ,a„] = ^ k(0^, vr) [oi, . . . , a„]. 



However the vanishing of mixed moments means now that the only tt 
which contribute are those which do not connect elements from different 
sets. Furthermore, the fact that all our variables are centered excludes 
singletons. But then it is easy to see that there are no such vr at all, so 
the sum is zero. 

2) Now we have to consider two cyclically alternating and centered 
tuples ai, . . . , and bi, . . . ,bn- Then we have 

(P2{ai ■ ■ - am] bi- ■ - bn) = ip{lm+n,1m.,n)[ai, • • • , ^m, &n] 



Again, the vanishing of mixed moments requires that (V, tt) connects 
only elements from the same set and the centredness of the elements 
excludes singletons. It is then easy to see that, for n > 2, the only 
possibilities for such (V, vr) arise for m = n and they have to be disc 
permutations (0,r, tt) which are pairings (ai, 6i+s)(a2, 62+s) ■ ■ ■ (on, ^n+s) 
for some s. The factors k{l2, {...))[ak,bk+s] are just <y?i(a262+s), so that 



7r)[ai, . . . ,a„] = 



TT£NC(n) 
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one finally gets, for n > 2, the formula 



n 



V52(ai ■ ■ - am] hi 



k=l 



'n+s 



For n = m = 1 one gets with 



ip2{ai; bi) = k2{ai, 02) + h,i{ai; bi) 



the conclusion that if2{ci'i', bi) has to vanish if ai and bi are from dif- 
ferent sets. Nothing is required if both are from the same set. We see 
that we get exactly the defining properties for second order freeness 



3) It would be nice to be able to reformulate in a similar way the 
definition of higher order freeness in terms of the (f instead of the 
cumulants. However, the situation with more than two circles is getting 
much more involved and we are not aware of such a reformulation for 
third and higher order freeness. 

As in the case of the first order freeness, one sees immediately that 
constants are free from everything. 

Proposition 7.8. Let {A,ip) be a HOPS. Then {1} is free of all orders 
from every subset X d A. 

Proof. We have to prove that 



unless n = 1. We will do this by induction on n. The case n = 2 is 
clear because 



from IMSQU. 



'^(In,7n(l),...,n(r))[l,a2, . . . 




«:(l2,(12))[l,a2] = <^i(l-a2)-^i(l)-<^i(a2) =0 



and 



/t(l2,(l)(2))[l,a2] = ¥^2(l;a2) = 0. 
In general, one has 




(V,7r)6-PSivc(n(l), ...,"('■)) 
(V.T)l<l(ln,7„(l),...,„(r))l 
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By induction hypothesis, in the later sum exactly terms of the form 
({1} U V, (1) U tt) with 

(V,7f)GP5jvc(^(l)-l,n(2),...,n(r)) 

contribute. In the case n(l) > 1 the sum over those yields 

<^(ln-l,7n(l)-l,n(2),...,n(r)[a2, ■ ■ ■ , On] ■ 

In this case, also 

</'(ln,7n(l),...,n(r))[l, ^2, ■ ■ ■ , On] = </?(ln-l, 7n(l)-l,n(2),...,n(r-)[a2, ■ ■ ■ , On], 

and thus «;(ln) 7n(i),...,n(r))[l, ^2, • • • , On] = 0. If, On the other side, 
n{l) — 1 (i.e., 1 is the only element on its circle), then we have to 
set 

P5;vc(O,n(2),...,n(r)) = 0, 

because then the first circle cannot be connected to the others if we 
ask 1 to be a cycle of its own. But this means that in this case 

'«(ln,7n(l),...,n(r))[l, ^2, ■ ■ ■ , On] = '/'(In, 7n(l),...,n(r-)) [1, 02, ■ ■ ■ , On] 

However, for n{l) — 1 and n > 1 we have 

</'(ln,7l,...,n(r-))[l,a2, . . . , a„] = 0. 

□ 

Note that our definition of freeness behaves clearly very nicely with 
respect to decompositions of our sets. For example, we have that 
Xi,X2, are free if and only if Xi and X2 U A3 are free and X2 and 
A's are free. Thus we can reduce the investigation of freeness to the 
understanding of freeness for the case of two sets. A characterization 
for this is given in the next theorem. 

Theorem 7.9. Let {A, f) be a higher order probability space and con- 
sider two subsets of Xi, X2 C A. Then the following are equivalent. 

(1) The sets Xi, X2 are free of all orders. 

(2) The sets X\ U {1}, X2 U {1} are free of all orders. 

(3) We have 

(/3(W,7)[aifoi,...,a„6„] 
= X] /«(V,7r)[ai, . . . ,a„] • <yi3(>V,(T)[6i, . . . ,6„] 

(V,7r)-(W,a)={W,7) 

for all n e N, all {U, 7) G VS{n) and all ai, . . . , a„ e A*! U {1}, 
e X2l^{l}. 
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(4) We have 

(/7(W,7)[ai6i, . . .,anbn] 
= X] <^(V,7r)[ai, ...,a„] •«;(>V,(7)[6i,...,6„] 

(V,7r)-(W,o-) = (W,7) 

for all n eN, all (U, 7) e VS{n) and all ai, ... ,an & XiU {1}, 

bi,...,bne X2U{1}. 

(5) We have 

K{U,-f)[aibi, . . .,anbn] 
= K(V,7r)[ai, ...,a„] •«;(W,(7)[6i,...,6„] 

(V,7r)-{W,o-) = (W,7) 

for all n eN, all {U, 7) e VS{n) and a// oi, . . . , a„ e A*! U {1}, 
6i,...,6„ G A'2U{1}. 

In order to prove this we would like to write ip{U, 7) [ai6i, . . . , a„6„] 
in the form <^(W, 7)[ai, 61, . . . , a„, Let us introduce the following 
formalism for this. Let (^,7) G VS{n) be a partitioned permutation 
of the numbers 1, 2, 3, . . . , n. Double now this set of numbers by intro- 
ducing a copy i, 2, 3, . . . , ri and interleave the new and old numbers as 
follows: 

1, i, 2, 2, 3, 3, ... , n, n. 

If we induce now (^^,7) on 1,2, ...,n to (^,7) on l,i,...,n, n by 
putting 

7(/c) = k and 7(^) = ^{k), 
then this has exactly the wanted property. The vanishing of mixed 
cumulants means that in the factorizations of (W, 7) in (V, tt) times 
a disc permutation we are only interested in (V, tt) which have the 
property that each block of V contains either only unbarred numbers 
or only bared numbers, i.e., (V, tt) must be of the form (VaUVb, TTaUTTb) 
with 

(V„, 7r„) e VS{1, ...,n) and (H, e . . . , n). 

Let us first observe some simple relations between the quantities on 
1, . . . ,n and their relatives on 1, 1, . . . , n, n. 

Lemma 7.10. 1) We have 

|7|=n+|7|, \l(\=n+\U\, 

and thus 

|(ZY,7)|=n+|(W,7)|. 
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2) We have 

|7r„ UTTfol = iTTal + \7lb\, \Va U Vb\ = \Va\ + \Vb\, 

and thus 

I {Va U Vb, 7r„ U TTfe) I = I (V„, TT.) I + I {Vb, TT,) | . 

3) We have that (vTa U7rfe)7 maps unbarred to bared and bared to unbarred 
elements and, for all k — 1, . . . ,n, 

[(TTa U7rfo)7]^(A;) = irbTValik), 

thus 

KTTaUTTb)-/! =n + |7rfe7r„7| 
Proof. Only the third part is non-trivial. To see this observe 

{TTa U 7rb)^{k) = TTainik)) 

and thus 

[(7raU7r5)7]'(^) = 7rfe7r„(7(A;)), 

which is our first equation, with the identification of tt^ e S{\, .... I) 
with the corresponding permutation in . . . , n). Since the mapping 
between bared and unbarred elements is clear, this yields that (tTo U 
■Kb)! and TTbiTaj have the same number of orbits which gives the last 
equation. □ 

This lemma allows us to characterize the contributing factorizations 
in {U, 7) in terms of special factorizations of {U, 7). 

Proposition 7.11. The statement 

{VaUVb,7raU7rb) eVSNciU,^) 

is equivalent to the statement 

{Va,Tra)-{Vb,Trb)eVSNc{U,j), 

where in the last product we identify (Vfe,7r) G VS{1, . . . ,n) with the 
corresponding element in VS{1, . . . ,n). 

Proof. Note that (Va U V;,, tTq U TTf,) G VSnc^!^^ l) is equivalent to 

(54) |(V„UH,7r„U7rb)| + |(7r„U7rfe)-H| = |(W,7)| 
and 

(55) W = (V„UH)V7. 

On the other hand, (Va,7ra) • {yb,T^b) £ VSNciJ^,!), nieans 
(V„,7r„) • (H,7r,) • (0 i,-i^,7r,-V„-S) = (^,7), 
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which is equivalent to 



(56) \iVa,na)\ + |(H,vr,)| + ^,-^-^1 = |(W,7)| 
and 

(57) W = V„ VVfe V 7- 

Equations and are, by Lemma 17. 101 equivalent. 

The equivalence between and (jSZ|) is also easily checked. □ 

Equipped with these tools we can now prove our main Theorem 17.91 

Proof. The equivalences between (3), (4), and (5) follow by convolving 
with the C or the fi function. That (2) is actually the same as (1) 
follows from Prop. 17.81 
(1) =^ (3): We have 

V9(W,7)[ai6i, . . . ,a„6„] = ip(U,^)[ai,bi, . . . ,a„,6„] 

= ^ K(V,7r)[ai,6i, . . . ,a„,6„] ■ C(>V,cr) 

(V,7r)-(W,o-) = (W,7) 

= ^ K{V,n)[ai,bi, . . . ,an,bn] 

(V,7r)e-P5]vc(W,7) 

By our assumption on the vanishing of mixed cumulants, only (V, vr) 
of the form (Vi U V2, tt^ U Hb) with 

(Va, Tfa) e VS{1, ...,n) and (H, vr^) G VS{1, ...,n) 

contribute and, by the above Proposition I7.1H 

is equivalent to 
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Thus we can continue with 
V?(W, 7) [ai6i, . . . ,a„6„] 

= ^ K{Va,7ra)[ai,a2, . . . ,an\ ■ K(yb,TTb)[bi,b2, . . . ,bn] 

(V„UV6,7r„U7rf,)eP5jvc(W,7) 

= ^ K{Va,7ra)[ai,a2,...,an\-KiVb,n)[bi,b2,...,bn] 

= ^ K{Va,7ra)[ai,a2, . . . ,an]- 

(Va,7ra)-(V6,7r6)-{W,a) = (W,7) 

■/€(H,7rfc)[6i,62,...,&n]-C(W,^) 
= ^ K{ya,7!-a)[ai, . . . ,an] ■ ip(y,n)[bi, . . . ,bn]. 

(Va,7ra)-(V,7r) = (W,7) 

(3) (1): Note that (3) allows us to calculate all moments of ele- 
ments from U X2 out of the moments of elements from Xi and the 
moments of elements from A2. (In order to do so, we also have to allow 
some of the a's or 6's to be equal to the unit 1.) Since this calculation 
rule is the same as for free sets, this shows that the sets Xi and X2 
must be free. □ 

This theorem is now the key ingredient to transfer freeness from sets 
to their generated algebras. 

Theorem 7.12. Let {A,(f) be a HOPS and consider subsets 

For each i & I , let Ai be the unital algebra generated by elements from 

Xi. Then the following are equivalent. 

(1) The subsets {Xi)i^i are free of all orders. 

(2) The subalgebras (^j)jg/ are free of all orders. 

Proof. Since the cumulant /t(V, vr) [ai, . . . multi-linear func- 

tional in the n variables it is clear that taking sums of el- 

ements within the sets Xi preserves freeness. What we have to see is 
that also taking products preserves freeness. Since we can iterate our 
arguments, it suffices to see the following: if Xi and X2 are free, then 
also Xi U {aotti \ oq, ai G Xi\ and X2 are free. Adding one product 
after the other to Xi and by Theorem 17.91 it is enough to show that 

V9(W,7)[aoai&i,a2&2 • • • 

= XI V5(V,7r)[aoai,a2 . . . ,a„] ■ k(W, (t)[6i, . . . ,6„] 

(V,7r)-{W,a) = (W,7) 

for all n G N, all (^,7) G VS{n) and all a^.ai, . . . ^ttn E XiU {1}, 
61, . . . , 6„ G A's U {1}. Let us induce (W, vr) G . . . , n) to (W, tt) G 
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VS{0, 1, . . . , n) by requiring that VV and vr restricted to 1, . . . ,n agree 
with W and vr, respectively, and that and 1 are in the same block of 
VV and 7r(0) = 1. Then we can calculate 

V?(W,7)[aoai&i,a2&2 • • • , ctn&n] = V5(W, 7r)[aol, ai&i, 02^2, • • • , fln&n] 
= ^ v5(V,7r)[ao,ai,a2 . . . ,a„] ■ fi:(W, (t)[1,6i, . . . 

(V,7r)-(>V,o-)=(W,7) 

By Proposition 17 . 81 we know that ft(>V, cr)[l, &i, . . . , bn] is only different 
from zero if W has as a singleton, i.e., (W, a) has to be of the form 

W = {0}UW, a={0)a, 

with 

GP5(l,...,n). 

But then we must have that 7r(0) = 1 and and 1 must be in the same 
block of V. Thus there is a unique (V, vr') so that (V, vr) = (V', vr') and 

(V,7r)-(>V,a) = (W,7) 

is equivalent to 

(V',7r')-(W,a) = (W,7). 
Note also that in this situation 

K(yV,cr)[ao,ai, , . . . , a„] = k(VV, 62, • • •, inl- 

and 

ipiy', vr') [ao, ai, , . . . , a„] = ip{y', 7i')[aoai, 02, ... , On]- 
So we can continue the above calculation as follows 

V9(W, 7) [aoai&i,a2&2 • • • 

= X] V5(V',7r')[aoai, a2 . . . ,a„] ■ ft(VV, a-)[6i, . . . 

(V',7r')-{W,S) = (W,7) 

which is exactly what we had to show. □ 

7.4. Distribution of one random variable. For the case where we 
restrict our attention to just one random variable a e .4 we introduce 
the following notation. 

Notation 7.13. Let (^, v?) be a HOPS and let a e A. 
1) For, (V, vr) G VS{n), we will write 

^'^{V,n) :=y.(V,7r)[a_^ 

n-times 
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and 

K^iV^ir) := /t(V, tt) [a, . . . ,a] . 

n-times 

2) A Young diagram is a A = (Ai,...,Ai) for some / G N and 

Ai, . . . , A/ e N with Ai > A2 > ■ ■ ■ > A;. We put |A| := Ai H h A/ 

(the total number of boxes of the Young diagram A). The set of all 
Young diagrams will be denoted by Y. 

3) The information about the higher order moments of a can also 
be parametrized by Young diagrams as follows: for A = (Ai, . . . , A/) we 
put 

V^'^(A) := ipilix\, vr) \a^_^^ = v^/(a^\ . ..,a^') 

n-times 

where vr is any permutation whose conjugacy class corresponds to A 
(i.e., vr G ^ia] has cycles of length Ai, . . . , A;. The collection of all higher 
order moments {f^{^))xeY called the (higher order) distribution of 
a. 

4) Similarly as for moments, we put 

k"(A) := ft:(l|A|,7r) Ja, . .^.,a] , 

n-times 

where vr is any permutation whose conjugacy class corresponds to A. 

Remark 7.14. For first and second order moments and cumulants, we 
used in Section |21 also the following notations: 

an := f"'{ln, In) Oi^^^n '■= V5"( Wn, 7m,n), 

and 

where 7^ and 7^,^ are permutations with one cycle and two cycles, 
respectively. 

The vanishing of mixed cumulants translates in this framework into 
the additivity of the cumulants for sums of free variables. 

Theorem 7.15. Let {A, (p) he a HOPS and a,h E A free of all orders. 
Then we have 

k''+\\) = k\\) + k\\) 



for all A G Y. 
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Proof. By the multilinearity of the cumulants and the vanishing of 
mixed cumulants for free variable, we have for any n G N and vr G 5*^,: 

K^+\ln,7T) = K{ln,7r)[a + b,...,a + b] 

= K{ln, Tc)[a, . . . , a] + K{ln, tt) [6, . . . , 6] 

□ 

8. Random matrices, Itzykson-Zuber integrals and higher 

order freeness 

8.1. Asymptotic higher order freeness of random matrices. Let 

us now come back to our original motivation for our theory - the as- 
ymptotic behavior of random matrices. In order to reformulate our 
calculations from Section 0] in our language of higher order freeness we 
still need to define the notion of "asymptotic freeness" . 

Definition 8.1. 1) Let {A,(f) and, for each G N, (An,^^'^'^) be 
HOPSs. Let I be an index set and for each i G /, Oj G ^ and a[^'' G An 
(A^ G N). We say that the family (a^-^'' | i G /) converges, for ^ oo, 
to (ttj I i G /), denoted by 

if we have for all n G N and all polynomials pi, . . . ,p„ in |/|-many 
non-commuting indeterminates that 

(58) limAT^oo V^i^^ (^pi((aj^^)ie/), . . ■ ,Pn{{a'l^^)iei)^ 

2) Let, for each G N, (^^v, ^^^^) be HOPSs. Let I be an index set 
and, for each i G / and A^ G N, a-^-* G An- We say that the sequence 
of families (a-^'*)jg/ has a limit distribution of all orders if there exists 
a HOPS {A, ip) such that 

[o-i )iei [o-ijiei^ 

for some tti E A {i E I) 

3) Let, for each A^ G N, {An,'^^^^) be HOPSs. Let / be an index 
set and, for each i G / and A^ G N, a-^-* G An- Let J = Ji U • • • U 
be a decomposition of / into k disjoint subsets. We say that the sets 
{a\'^^ I i G Ji}, . . . , {a\'^^ \ i G Ik} are asymptotically free of all orders 
if there exists a HOPS {A, ip) such that 

r W\ ( \ 

[O-i )i&I {O'iji&I-i 
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for some ai E A {i E I) and such that the sets {oj | i G Ji },..., {oi | 
i e Ik} are free of all orders in {A, ip). 

With this notation and by invoking Theorem 17. 91 we can reformulate 
our main result on random matrices, Theorem 14 .41 in the following 
form. 

Theorem 8.2. Let Mn ■= Mn ® L°°-{Q) be an ensemble of N x N- 
random matrices. Define rescaled correlation functions ip^^^ = (<^n^'')neN 
on Mn by (n eN, Di, . . . ,Dn e Mn) 

(59) ^L'^HDi, • • • , ^n) := A;„(Tr(Di), . . . , Tr(D„)) ■ iV^"". 

Assume that we have, for each N eN, subalgebras An, ^n ^ Mn such 
that 

(1) An is a unitarily invariant ensemble, 

(2) An and Bn are independent. 

Let be a family of elements in {An-, <f^^^) which has a higher 

order limit distribution and let {Bj^^)j(zj (N G be a family of ele- 
ments in {Bn,'^^^^) which has a higher order limit distribution. Then 
the families {A-^"* | i G /} and {-B]^'* | j G J} are asymptotically free 
of all orders. 

8.2. Itzykson-Zuber integrals. 

Definition 8.3. For N x N matrices An, Bn their Itzykson-Zuber in- 
tegral is defined as the following function in 2; G C: 

IZ{z,An,Bn) := Ar-2logE(e"^^(^^^^^^*)), 

where U denotes a Haar unitary N x A^-random matrix. 

Consider now a sequence of such matrices An and Bn- Note that 
An and Bn are non-random, thus all distributions of order higher 
than 1 vanish identically. If we assume that An and Bn have a first 
order (eigenvalue) limit distribution for N 00, then it is known (see 
|Coin8j ) that each Taylor coefficient about zero of z — >• IZ{z, An, Bn) 
admits a limit as — > 00. Note that the effect of the Haar unitary 
random matrix in the above Itzykson-Zuber integral was to make An 
and UBnU* asymptotically free of all orders. We show now that this 
kind of result extends also to the case of random matrices An and Bn, 
and that our theory allows to identify the limit of the Taylor coefficients 
very precisely. 

Theorem 8.4. Let A = {AN)NeN and B = {Bn)n<^n be two ensem- 
bles of N X N -random matrices which are asymptotically free of all 
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orders with respect to the rescaled correlation functions (p^^\ Denote 
the corresponding limiting distribution of {Aj^)]\r^^ by and the cor- 
responding limit distribution of {Bfyf)j^^fq by if^ . Then, as formal power 
series in z, we have 
(60) 

oo „ 

lim Ar-2logE[e^^^(^^^~)] ^ /s:'^(V, vr) -/(W, a). 

n=l ^' (V,7r),(W,CT)6PS(n) 
(V,7r).{W,<T) = (l„,e) 

Proof. Recall that the logarithm of the exponential generating series of 
the moments of a random variable is the exponential generating series 
of the classical cumulants of that variable. Thus we have 

. logE[e"^^^^'^] 



z 

V fc„(iVTr(Ajv5^), . . . , NTi{AnBn)) ■ -r 
^ — ^ n\ 

n=l 



n=l 

By our assumption that Ajq and -Bat are asymptotically free with re- 
spect to (pi^'^ = iV"~^(y9^, this converges to 

°° 

Vy?(l„,e)[a6, ...,a6] ■ — , 
^-^ n\ 

n=l 

where a and b are free of all orders with respect to ip. Theorem 17.91 
yields then the assertion. □ 

The main result of this part is the following theorem, which shows 
that the higher order analogue of the Itzykson-Zuber integral behaves 
like a kind of i?-transform in one matrix argument if the other is re- 
stricted to non-random matrices. This result can be seen as a higher 
order version of a result of Zinn- Justin |Z.T99j . 

Theorem 8.5. Let C = {Cn)^^^ be a sequence of non-random N x 
N matrices which has a first order limit distribution. Then, for any 
sequence of N x N -random matrices A = {An)^^^ for which a higher 
order limit distribution exists, we define as a formal power series in the 
moments of C 

R^{C):= lim iV-2logE[e^^(^^^^^^^^)], 

where Un are Haar unitary N x N random matrices which are inde- 
pendent from A]\f. Then we have the following: 
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1) If A = (Aiv)ArGN ond B = (i?Ar)7VGN 0,1"^ random matrix ensembles 
which are asymptotically free, then we have as formal power series 

(61) /?^+^(C) = R'^iC) + i?^(C). 

2) More precisely, if we denote the limit moments of C by Xk '■= 
liniAr^oo ti'(C*Ar) (k eN), then one has as a formal power series 

(62) i?^((a;,),^^) = ^^^A^l|,|,A), 

where c\ is the number of permutations in the conjugation class in S\x\ 
corresponding to X, and where 

x^ ■=x^'x2^-- - , for \={\i>\2 >...). 

Proof. In order to get ()62|) (which imphes ()6H)) we have to speciahze 
(I60|) to the situation that B = C has all moments of order higher than 
1 equal to zero. This means that ^p^iW, cr) = unless W = Oo-, in 
which case we have 

if\Q„,a)=x\ 

where A is the Young diagram corresponding to cr. (jH^ follows then 
from the simple observation that, for fixed cr, the only solution of (V, vr)- 
(0^, cr) = (1„, e) is given by (V, vr) = (1„, a"^). □ 

This theorem tells us that Itzykson-Zuber type of integrals contain 
the whole data to linearize higher order freeness. Also, the technology 
introduced in this paper, especially Theorem 16.31 gives methods to re- 
fine asymptotics of spherical integrals a la Guionnet and Maida (see 

IGMOni). 

As a foreshadowing of such applications we close with the following 
proposition. More details will be provided in a forthcoming paper. 

Let Cat be a non-random matrix of rank 2 with eigenvalues x and y 
(which will be considered as indeterminates in the following). Then we 
have that ip^^^ (W, a) [C, . . . , C] is only different from zero for W = Oo-, 
in which case it is 

^(^)(0., cr) [C, . . . , C] = (x^^ + 2/"^) ■ ■ ■ (x^' + y^O, 

where A(cr) := A = (Ai,...,Ai) is the Young diagram encoding the 
conjugacy class of a. Thus Theorem 14.41 vields in this case 

iV-MogE[e™-^^-^*] = ^Ar"-V(^)(ln,e)[A^C^,...,AA.C^] ■ ^ 

n=l 

oo n 

= E ^ E '^^'"Hin, <y-')[AN, . . . , A^] ■ {x'^^"^ + y'^^''^) ■ ■ ■ ix'^^"^ + y''^"^) 

71=1 ' a£S„ 
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Now we invoke the assumption on the existence of a hmit distribution 
for A]\f and the fact (which we have never used up to now) that the 
orders in the Weingarten function and thus also in k^^^ decreases in 
steps of 2. This allows us to recognize the two leading orders of the 
above quantity and we get the following result. 

Proposition 8.6. Let Cn be a matrix of rank 2 and eigenvalues x,y. 
Consider a sequence of N x N -random matrices A = (y47v)jveN for 
which a higher order limit distribution ip"" exists. Denote the limiting 
first and second order cumulants by 

where 7„ and 'ym,n ore permutations with one cycle and two cycles, 
respectively. Then, as a formal power series in x, y we have 

hm logE[e^^^-^^-^*] 

n ^-^ mn 

n>l m,n>l 

This expansion extends results of [ Col03j . It is more general than 
that obtained in jUMOSj, because it also handles the case when Ajsr has 
asymptotic fluctuations. However, unlike in |GM05j . the convergence 
that we obtain is formal and it would be very interesting to check if it 
still holds at an analytic level. 

9. Appendix: Surfaced permutations 

In this appendix we will present a more geometrical view on par- 
titioned permutations. As we shall see in the following, partitioned 
permutations are just special cases of "surfaced permutations"; in par- 
ticular the results of this article can be equivalently formulated in the 
language of surfaced permutations. On the other hand, for the pur- 
pose of this article we do not need anything more than just partitioned 
permutations and the Reader not interested in surfaced permutations 
may skip this Section without much harm. 

9.1. Motivations. Our goal is to study factorizations of permutations, 
i.e. solutions (tti, . . . , tt^) of the equation 

7 = TTl ■ ■ ■ TTfc, 

where 7 G S'„ is some fixed permutation and tti, . . . , tt^ G Sn are sub- 
ject to some additional constraints, depending on a particular context. 
Typically, one of these constrains concerns |7ri| + ■ ■ ■ + |7rfc|, the other 
one concerns the orbits of the action of tti, . . . , vrfc. 
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Figure 1. Example of a surfaced permutation. Its 
support is equal to (1,3)(2)(4) G S^. This surfaced 
permutation corresponds to a partitioned permutation 
({1,3,4}{2},(1,3)(2)(4)) 

It would be very useful to equip permutations tti, . . . , vr^ with some 
additional structure in such a way that the product tt i ■ • • vr^. of the 
resulting enriched permutations vfi, . . . , vr^ would carry both the infor- 
mation about the product vti ■ ■ ■ vTfc of permutations and the informa- 
tion about IttiI + ■ ■ ■ + |7rfe|. As we shall see in the following, surfaced 
permutations provide an appropriate tool. 

9.2. Definition. We say that a = {S,j) is a surfaced permutation 
of some finite set A if S" is a two-dimensional surface with a fixed 
orientation and with a boundary dS and if j : A ^ dS is a injection. 
We can think about the information carried by j as follows: some of the 
points on the boundary dS are distinguished and carry different labels 
from the set A. We also require that every connected component of 
dS carries at least one distinguished point. An example of a surfaced 
permutation is presented on Figure [T] 

We identify surfaced permutations (S'i,ji), (5*2, j2) of the same set 
A if there exists a orientation preserving homeomorphism / : 5*1 — > 5*2 
such that foji = j2. The set of surfaced permutations of set {1, . . . ,n} 
will be denoted by SSn- 

9.3. Surfaced permutations and the usual permutations. Let 

{S,j) G SSn', the boundary dS with the inherited orientation from S 
is just a collection of oriented circles with some distinguished points 
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Figure 2. Convention for splitting labels. 

labeled 1, . . . , n marked on them. In this way we can define a permuta- 
tion a G S'„, called the support of {S,j), the cycles of which correspond 
to connected components of dS, as it can be seen on Figure H It is 
therefore a good idea to think that a surfaced permutation is just a 
(usual) permutation a E Sn equipped with some additional informa- 
tion carried by the surface S. 

A surfaced permutation {S,j) G SSn can be uniquely specified (up 
to the equivalence relation) by its support a E Sn and by specifying 
the shape of the connected components of S. The latter information is 
given by an equivalence relation on cycles of a (each class corresponds 
to a connected component of S) and furthermore for each class of this 
relation we should specify the genus of the corresponding connected 
component of S. Above it should be understood that the genus of 
a surface S with a boundary is by definition equal to the genus of a 
surface 5" without boundary obtained from S by gluing a disc to every 
connected component of dS; for example both a disc and the lateral 
surface of a cylinder have genus zero. 

9.4. Surfaced permutations and partitioned permutations. 

Among surfaced permutations a special class will be very important 
for our purposes, namely surfaced permutations {S,j) such that each 
connected component of S has genus zero. It is easy to see that there is 
a bijection between such surfaced permutations {S,j) and partitioned 
permutations (V, a) given as follows: a is the support of {S,j) and V 
is the partition given by connected components of S. 

9.5. Products of surfaced permutations. Let surfaced permuta- 
tions (S*!, ji), (5*2, ^2) G <SSn be given. On the boundary of 5*2 there are 
marked points labeled by numbers 1, . . . ,n; let us spht every marked 
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point k into a consecutive pair of points k and k', as it is presented on 
the example from Figure |21 In the second step, for each k G {1, . . . ,n} 
we glue a small neighborhood of the vertex k G dSi to a small neigh- 
borhood of the vertex k' G 882 in such a way that the orientations of 
5*1 and 5*2 coincide. In this way we obtain a new surface S which has 
marked points on its boundary dS and these are exactly the vertices 
from 882 labeled 1, . . . ,n; we denote the resulting surfaced permuta- 
tion by {S,j) and we call it a product ("S"!, ji)(S'2, ^2) of the original 
surfaced permutations. This choice of gluing surfaces Si and 5*2 im- 
plies that the support of ji)(S'2, ^2) is equal to the product of the 
support of (5*1, ji) and the support of (5'2, j2)- 

It is not difficult to explain now the definition of the product of 
partitioned permutations (Definition 14. 9|) : we treat partitioned per- 
mutations as surfaced permutations and compute their product; if the 
genus of the resulting surface is zero we can identify it with another 
partitioned permutation, otherwise we set the product to be zero. 

It is not difficult to show that for surfaced permutations the prod- 
uct is associative and the associativity of the product of partitioned 
permutations is a simple corollary. 
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